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RMC, PMC, YMC 
^ ^ 

DjYj L ^ 

RM,PM,YM 

DU,YU,LU, 

RU,PU,YU 
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fc,f 

ibJ b ikb 


projected frontal area of container (eq. 4.6) 

corrections of YC, LC and YMC to theoretical zero values on 
boundaries of R 0 

these denote the wind axes components of the aerodynamic force and 
moment vectors, which are (— D,Y,—L ) and ( RM,PM,YM ) , 
respectively: termed drag, sideforce, lift, and roll, pitch, yaw moments, 
respectively. Moments are measured about the geometric center of the 
box. 

the NASA Ames Research Center’s wind tunnel measurements of 
D,Y,...,YM (ref. 3) corrected for interference 

measurement errors, DC — D,YC — Y, etc. 


estimated values of the force and moment components obtained in 
this paper 

wind tunnel measurements from University of Maryland (ref. 6) 


attitude domain of interest in simulation: —180° < ip < 180°, 
—90° < a < 90°, and its quadrants (fig 3.1). Do is the 
quadrant 0 < rp, a < 90° containing the measurements. 

elementary rotation matrices for rotation through angle a about 
i,j_ and k axes, respectively. 

{E 1 (a),E 2 (a),E 3 (a)} = 

0 0 \ / cos o 0 — sincr\ ( cos a 

cos o sin o I , I 0 1 0 I , I —sin a 

— sin a cos a J ^sino 0 cos a J \ 0 

two-dimensional vectors, ( YC,LC) T and (Y, L) T 

container body axes; orthogonal right-handed frame with i b 
parallel to a long edge and perpendicular to the bottom face 
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Sl,S 2 , S 3 , 
S 4 , S 5 , S& 


tunnel mount axes defined from the angle measurements (ip, a), 
where 26 0(U , 0 = E 2 (a)-&3(0) 

wind axes; orthogonal right-handed frame with t w along Va. and k w 
in the body vertical plane. 

tunnel axes in which forces and moments are measured; orthogonal 
right-handed frame with i w and k Wo along the tunnel longitudinal axis 
and the local vertical, respectively. 

dynamic pressure, pVa 2 / 2 

half of the sum and the difference, respectively, of the roll moment 
measurements at pairs of equivalent points (eq. 4.14) 

regions in V 0 bounded by the lines a = 45°, 

{V>*(a), ip( 45 )(a), 0 < a < 90°}, respectively (figures 3.5 and 3.6); 

R. 0 is the region spanned by the measurements. 

the subdivision of V 0 defined by the lines a = 45°, {ip* (a)} , {ip^ 45 ^ (a)} 
(fig. 3.6). 


Sl,S2,S3,S 4 , 
S5, Sq, S7 


^a,6 


Vo, Va 


diagonal matrices which define symmetry properties and correspond to 
rotations of 180° about various axes (fig. 5.1). 

transformation of a vector from reference frame “b” to frame “a” , 
where the subscripts a and b can be any of b,b 0 ,w, or w 0 to indicate 
the reference frames used in this work (eqs. 3.1 - 3.3 ) 

air-velocity vector, airspeed 


a angle of attack; tan a = Va • k^/Va • t 6 , —90° < a < 90° 

P sideslip angle; sin/? = Va, • j /V o, sign(cos/3) — sign(Va • t 6 ); 

-180° <0< 180° 


6D,6Y,SL differences between values of D, Y, and L given from the Ames 

Research Center’s data at equivalent points 

AD, AY, A L, interference tares (appendix A) or the differences between data from 

ARM, A PM, AYM references 3 and 6 (appendix B) 
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angle between air velocity vector and longitudinal body axis (eq. 4.7) 
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<Pw 
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Via) 

0< 45 )(a) 


0(0, a) 
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()T 


roll angle locating wind axes relative to tunnel axes. T w>Wo = E i(4> w ) 

model heading angle relative to tunnel or wind axes: 0 = —0 

heading angle along the line of self-equivalent points (fig. 3.6) 

heading angle along the line equivalent to the boundary of £ 0 , 
a = 45° (fig. 3.6) 

heading angle such that £(0,45°) = £(0, a) (eq. 4.8) 
vectors are indicated by an underbar 

error in the estimate or measurement of ( ) 
estimate of ( ) 
transpose of ( ) 


UNITS 

English units (lb, ft, and sec) and degrees are used throughout and are generally omitted in the 
figures and tables. Aerodynamic forces and moments are measured as force or moment divided 
by dynamic pressure, which have units of ft 2 and ft 3 , respectively. 



SUMMARY 


A comprehensive model of the static aerodynamics of the 8- by 8- by 20-ft container known 
as the MILVAN is determined by combining the available wind tunnel data from a 1972 NASA 
Ames Research Center study taken over the restricted attitude domain, {0 < ip < 90°; 

0 < a < 45°}, with an extrapolation to the complete domain of interest in simulations, 

{ — 180° < xp < 180°; —90° < a < 90°}, based on an appropriate theory derived from the 
geometric symmetry of rectangular boxes. The theory predicts the symmetry properties of all 
force and moment components about various boundaries, including all zero lines, their behavior 
along certain boundaries, and their interrelationships at equivalent points. It was found that the 
aerodynamics at any attitude can be defined from the aerodynamics at an equivalent attitude 
in the restricted domain {0 < xp < 45°; 0 < a < 90°}. However, a similar comprehensive 
equivalence with the domain spanned by the measurements is not available; in particular, about 
two-thirds of the domain with |a| > 45° is unrelated to the domain of measurements. Nevertheless, 
an estimate can be given for this region consistent with the measured and theoretical values 
and the theoretical equivalence relations. Discrepancies between the data and the predicted 
symmetry properties ranged from small to gross, depending on the component, and were removed 
to apply the extrapolation theory. Small discrepancies were treated by imposing the required 
zero lines and averaging data at equivalent points. The roll- and pitching-moment data showed 
gross discrepancies and were discarded. Roll moment could be estimated using the available data 
from an independent study or assumed negligible, while pitching moment was estimated using the 
yaw moment data and the theory. The tunnel errors or unrealistic theoretical assumptions which 
account for these discrepancies cannot be formally identified because records of the original wind 
tunnel sensor readings and data processing had been discarded prior to the present study, nor were 
new tests within the scope of this study. However, limited data from independent wind tunnel 
studies are reviewed and show good to fair agreement with both the theory and the estimate given 
here. 


1. INTRODUCTION 

In rotorcraft slung load operations, a large variety of objects are carried which have a corre- 
spondingly large range of aerodynamic characteristics, and equilibrium and dynamic behavior in 
flight. The 8- by 8- by 20-ft container, also known as the MILVAN, is a standard cargo carrier 
and is a common load in these operations. It develops significant aerodynamic specific force and 
moment at low speeds, including static, quasi-static, and unsteady flow terms (refs. 1-8), and 
can become unstable at speeds well below the power-limited speed of the helicopter-load system. 
In addition, its simple shape results in static aerodynamics which are more easily documented in 
the wind tunnel and represented in simulations than those of arbitrarily shaped loads. Thus, it 
has served as the typical difficult load in much of the slung-load stabilization research of the past, 
and is the usual load in slung-load simulations (e.g., refs. 1 and 2). 

Various slung-load configurations are illustrated in figure 1.1, including the single point sus- 
pension, the two-point suspension for tandem-rotor helicopters, and the dual-lift system. The 
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MILVAN angle of attack in steady flight can range from —90° at very low speeds to small positive 
values. For the single point suspension, the MILVAN trims broadside to the flow at low speeds, 
and becomes unstable at 40 - 60 knots. In the two-point suspension the MILVAN trims lengthwise 
to the flow in its minimum drag orientation and is stable to much higher speeds. 

Future slung-load research may focus on the development of automatic stabilization and 
flightpath control systems, and the control of exotic suspensions such as the dual-lift system. 
These efforts can benefit from a comprehensive simulation model which captures realistic load 
motion for both static and dynamic condition. Towards this end, reference 8 provides a recent 
derivation of the system equations of motion for single-point suspensions, and the load kinematics 
in a two-point suspension are derived in reference 9. The present work adds a comprehensive 
model of the MILVAN’s static aerodynamics derived from the available wind tunnel data. This 
suffices for realistic calculation of the MILVAN’s attitude and the suspension forces in equilib- 
rium. However, as demonstrated in references 4 and 7 , the MILVAN’s dynamic behavior depends 
on significant unsteady aerodynamic effects, and a linear perturbation model of its unsteady 
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aerodynamics is given in reference 8 from the available wind tunnel. 

Wind tunnel measurements of the MILVAN’s static aerodynamics are available from several 
sources. The most comprehensive measurements were made at NASA Ames Research Center 
(ARC) and are documented in reference 3. These data (termed the Ames data hereafter) are 
the basis of the present work, but are almost entirely restricted to angle of attack between 0 and 
45° and heading angle between 0 and 90°. Additional independent data are available from an 
Army/Northrop study (refs. 4 and 5), and a VERTOL/ University of Maryland study (ref. 6) 
(UOM data). The data in references 4 and 5 survive only in graphical form, are restricted to a 
few angles of attack (0, 10°, 28°), and not all components are given. The data in reference 6 spans 
the same domain as that in reference 3, but no data were taken at heading angles between 50° 
and 90° where the side force and yaw moment exhibit peak magnitudes and strong nonlinearities. 
The reference 6 data were the source of the abbreviated MILVAN aerodynamic model in the ARC 
simulation of reference 2. Both of these studies provided important supplementary information 
in the present work. In addition, an Army /University of Bristol study, (ref. 7) provides results 
from surface pressure distribution measurements and flow visualization tests for a limited set of 
attitudes. 

The static aerodynamic model of the MILVAN in the Langley Research Center’s simulation 
(ref. 1) was also derived from the data used here. In that model, drag, lift, sideforce and yaw 
moment were given by trigonometric functions fitted to the data. This treatment is succinct, 
but of modest accuracy relative to the data in some regions because of the limited number of 
fitting function parameters, and erroneous in its extrapolation to angles of attack above 45° for 
lack of an extrapolation theory. The present work attempts to improve the static aerodynamic 
simulation model by; 1) providing tabulated function models, 2) deriving and applying a theory 
for extrapolating the model beyond the region of the measurements, and 3) adding the missing 
roll and pitching-moment components. Tabulated models can be used routinely and efficiently 
in digital simulations to represent arbitrary nonlinear functions to the accuracy of the available 
measurements . 

A theory for extrapolating the model beyond the region of measurements is derived from 
the geometric symmetry of the box. The box is indistiguishable to the flow and yields identical 
wind tunnel measurements when rotated 90° or 180° about its longitudinal axis or turned end 
for end, and also possesses the usual force and moment symmetries of bodies with horizontal 
and vertical planes of symmetry. These factors permit the aerodynamics at any attitude to be 
defined from the aerodynamics at a corresponding point in a small region. The smallest such 
region identified in the analysis is contained in the half-quadrant {0 < < 45°; 0 < a < 90°}. 

It is found that the half-quadrant spanned by the measurements of reference 3 does not suffice to 
define the aerodynamics at all attitudes. In particular, the aerodynamics in most of the region 
with angle of attack above 45° cannot be obtained from the measurements, but an estimate can 
be given consistent with the measurements along one boundary of this region, theoretical values 
along the remaining boundaries, and a theoretical relation between the aerodynamic components 
at equivalent points within the region. 

When the measurements of reference 3 are compared with the theory, discrepancies are 
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found which range from small for some components to gross disagreement in the case of roll 
and pitching moments. These discrepancies must be removed to apply the extrapolation theory. 
Small discrepancies are treated by imposing the theoretical zero lines and averaging the data at 
equivalent points. The roll and pitching moment data cannot be treated in this way because the 
discrepancies dominate the measured values and a corrected table would be completely uncertain. 
These data are discarded and roll moment can then be estimated from the available independent 
wind tunnel data or neglected, while pitching moment is estimated from the yaw moment data and 
the theory. The independent wind tunnel data from references 4-6 are also reviewed and show 
good to fair agreement with both the theory and the estimate derived from the Ames data for all 
components. The final model consists of data tables for all six force and moment components in 
wind axes and extrapolation formulas for both wind axes and body axes. 

We acknowledge the helpfulness of the authors and organizations of the wind tunnel studies 
(refs. 3 - 6) in locating and discussing the original data, and we regard with mixed feelings 
the penetrating questions of one reviewer which exposed an error in the fundamental symmetry 
equation in the first draft of this report that required all text and computations after Section 2 to 
be redone. Last, we note that all numerical analyses discussed here were carried out with a digital 
computer and that the tables of values in the text were printed directly from the computer. 
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2. A REFERENCE TABLE OF WIND TUNNEL DATA 

The Ames data, reference 3, consists of two tables of force and moment measurements uncor- 
rected for model support interference and two tables of interference tare measurements. The pairs 
of tables were generated by varying only angle of attack or only yaw angle in the sequences of 
measurements. The basic measurements were uncorrected for interference in reference 3 because 
of the erratic nature of the tare measurements which left their use open to individual interpreta- 
tion. These four tables are combined here into a single reference table of corrected measurements 
covering the region 

*o = {(0,«): 0€[O o ,9O°], a 6 [0°, 45°]} (2.1) 

The method and results are summarized next and details are documented in appendix A. 

First, the two tables of force and moment measurements are combined as follows. These 
tables cover the two regions shown in figure 2.1(a); these regions are not coincident but collectively 
they cover R 0 and can be averaged in the common region to reduce the standard deviation of 
random measurement errors there. The present treatment uses all available data in R 0 . An 
examination of the differences between the two tables shows; 1) small differences for the force 
and yaw moment components (below 10Q lb, 30Q ft-lb), which indicates good repeatability and 
small random measurement errors for these components, and 2) large differences for the roll and 
pitching moment data (25Q to 100Q ft-lb over most of R 0 ), which indicates large random errors 
for these components. 

Strong vibrations of the model and corresponding variations in scale readings were observed 
in the ARC tunnel tests, and the static aerodynamic data were obtained by averaging numerous 
readings. This average is expected to converge to a unique value if the number of readings is 
increased indefinitely, and to show only small differences among averages obtained from sufficiently 
large sets of readings. This is the case for the force and yaw moment data. The source of the 
much larger differences in the EM and PM data is unknown and cannot be determined without 
a more careful analysis of the test conditions and data processing than is now possible. In any 
case, the limiting average values of EM and PM have not been approached in the Ames data, 
and this is confirmed by the large differences between these data and the data from independent 
studies reviewed in appendices B and C, as well as from the values predicted by the theory of the 
next section, and by the general agreement among these other sources. 

Second, the support interference is estimated from the two tables of tare measurements as 
follows. The regions covered by the two tables are shown in figure 2.1(b); these are insufficient 
to cover R a . Further, the pitch sweep data measured the interference of only part of the support 
and are not used here except to indicate interference trends at angles of attack above 20° where 
there are no yaw sweep data. The yaw sweep data were reviewed with the object of estimating 
the interference function from the discernible systematic trends in these data and attributing the 
random variation around these trends to measurement errors. 

For the force- and yaw-moment components, the tare data showed consistent variations with 
and random variation with a; therefore, the interference was estimated by averaging the tare 
data over a. The averaged tares are small in magnitude, and the differences between the estimated 
tare and the original data are also small and statistically similar to the differences between the two 
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sets of measurements just discussed, as is expected if the same random error processes account for 
both sets of differences. The roll-moment tare data appeared to consist almost entirely of samples 
of a random process; the systematic interference was estimated as the table mean, which was 
nearly zero, while the data measured large interference moments which were statistically similar 
to the differences between the two sets of roll moment measurements just discussed. Evidently, the 
systematic roll moment interference was much smaller than the random-roll-moment measurement 
errors and not detectable in the Ames data. Finally, the pitching moment data showed discernible 
trends in both ip and a and the interference was estimated by ad hoc smoothing of the table of 
measurements. Results showed large estimated interference moments and smaller, moderate-sized 
differences between the estimate and the data. 

The resulting reference table is given in figure 2.2 and plotted in figures 2.3 and 2.4. These 
corrected data are indicated by appending “C” to the notation (DC, Y C, etc.) in the remainder 
of this paper. A cursory examination of the plots shows that drag increases with both ip and 
a, (with projected frontal axea) as expected, and that side force and yaw moment approximately 
pass through zero at ip = 0 and 90° consistent with their expected antisymmetry in ip about 
these points. However, lift shows significant nonzero values at a = 0 where antisymmetry in a is 
expected, and pitching moment shows a gross departure from the same expected antisymmetry. 
These discrepancies suggest the presence of systematic errors remaining in the data, which are 
small to gross depending on the component, and they are examined in greater detail in the 
application section of this paper. 

The remainder of the problem is to define the static aerodynamics at all possible attitudes 
from these data given on Z 0 - The extent to which this can be done based on the geometric 
symmetry of the box is examined in the next section, and the resulting extrapolation theory 
and the data are combined in the following section to obtain a comprehensive estimate of the 
MILVAN’S static aerodynamics. 


6 



60 


X 

□ 


PITCH SWEEP DATA 
YAW SWEEP DATA 


30 


□□□□□□□□ 

□□□□□□□□ 


xx>oo<xxxxx 

X 

X 

X 

XX 

XXXXXXXXXX 

X 

X 

X 

XX 

xxx>o<xxxxx 

X 

X 

X 

XX 

DBO 

B 

B 

BB 

xxxxxxxxxx 

X 

X 

X 

XX 

xxxxxxxxxx 

X 

X 

X 

XX 


BBGtSG H 



-30 i I I I I I I I 

(a) FORCE AND MOMENT MEASUREMENTS 


60 


X PITCH SWEEP DATA 
□ YAW SWEEP DATA 


30 


a 


X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 

X 


pna-Mnnwnnixn 

G 

Q 

G 

X X X X 




BPQQflPGgQQBQ 

G 

G 

G 

mrnBnDBGGBD 

G 

G 

G 

hjgghggbggbo 

G 

G 

G 

BEGGED Q HQ OEJQ 

G 

G 

G 


-30 1 I I I I I I I 

-90 -60 -30 0 30 60 90 120 

* 

(b) INTERFERENCE TARE MEASUREMENTS 


Figure 2.1.- Measurement grids of the Ames data. 


7 


a 

0. 

57.9 

60.3 

68.2 

77.5 

88.3 

96.8 

DC/Q 
109.6 121.6 

132.4 

140.9 

150.8 

165.5 

195.1 

213.4 

216.4 

2. 

60.9 

62.6 

70.3 

78.7 

89.4 

98.3 

111.4 

123.4 

134.8 

143.0 

152.1 

167.7 

195.9 

212.8 

216.4 

4. 

63.8 

66.1 

73.5 

80.0 

91.0 

101.1 

112.9 

125.3 

136.6 

144.8 

153.6 

170.3 

197.5 

213.1 

217.9 

6. 

66.8 

69.8 

77.2 

83.3 

93.2 

103.3 

116.3 

127.5 

139.0 

146.8 

156.3 

173.4 

198.0 

212.3 

217.6 

8. 

70.5 

73.1 

80.3 

87.1 

96.1 

106.0 

119.3 

130.2 

141.8 

149.9 

161.0 

177.6 

198.4 

212.1 

216.6 

10. 

74.2 

76.4 

83.8 

90.7 

100.4 

110.6 

124.2 

133.8 

145.1 

153.7 

165.0 

180.2 

199.1 

212.5 

216.9 

12. 

78.1 

80.0 

87.1 

94.4 

104.7 

115.9 

129.1 

138.0 

149.4 

157.6 

168.1 

182.1 

199.8 

211.5 

216.8 

16. 

85.7 

87.0 

94.2 

102.2 

112.2 

124.8 

139.1 

147.3 

156.1 

163.8 

172.9 

184.8 

200.8 

210.0 

216.5 

20. 

94.6 

97.1 

103.8 

111.3 

122.3 

133.6 

146.3 

152.4 

162.4 

169.1 

176.8 

187.3 

200.9 

208.5 

216.9 

25. 

107.0 

109.1 

116.4 

123.9 

133.3 

144.4 

155.3 

161.5 

171.0 

177.1 

184.2 

192.4 

202.1 

206.5 

216.2 

30. 

118.9 

120.4 

127.1 

136.1 

142.1 

152.6 

162.4 

170.5 

179.6 

185.2 

191.3 

197.2 

202.9 

206.0 

216.9 

35. 

126.6 

131.4 

139.1 

144.8 

149.7 

160.9 

170.5 

182.4 

189.6 

193.4 

198.7 

201.6 

204.6 

204.1 

216.1 

40. 

141.8 

143.3 

149.8 

154.6 

159.9 

166.7 

179.4 

190.2 

197.8 

202.4 

208.5 

208.6 

206.1 

204.2 

214.4 

45. 

157.2 

159.8 

163.8 

164.0 

166.5 

176.6 

188.9 

197.6 

205.9 

210.2 

215.9 

214.7 

209.6 

201.1 

215.2 

a 

0. 

-2.5 

10.1 

18.5 

24.3 

29.1 

38.4 

YC/Q 
53.5 63.2 

68.5 

73.7 

76.6 

80.0 

87.0 

70.6 

6.0 

2. 

-1.5 

11.2 

19.3 

23.9 

28.1 

37.5 

52.5 

60.8 

66.5 

70.5 

73.3 

78.0 

87.9 

70.7 

6.4 

4. 

-1.0 

11.6 

19.9 

27.8 

28.1 

37.9 

51.7 

60.3 

65.4 

68.1 

70.5 

76.2 

88.4 

71.6 

8.0 

6. 

-1.5 

12.4 

22.5 

28.0 

30.1 

40.0 

50.7 

58.7 

63.4 

66.8 

70.8 

77.8 

89.5 

72.4 

7.3 

8. 

-2.3 

12.8 

24.4 

28.8 

31.5 

41.5 

50.8 

57.3 

63.2 

66.3 

70.6 

78.7 

89.5 

73.4 

7.5 

10. 

-2.9 

12.8 

27.2 

31.3 

35.8 

44.5 

51.6 

57.3 

61.7 

65.8 

70.9 

78.6 

89.7 

74.2 

7.8 

12. 

-4.3 

13.2 

28.8 

33.2 

38.9 

46.2 

51.9 

56.3 

61.9 

66.0 

70.8 

78.6 

90.8 

75.1 

7.5 

16. 

-5.8 

12.9 

32.3 

42.1 

49.7 

51.4 

49.2 

53.5 

58.4 

64.2 

70.0 

77.1 

89.8 

76.9 

8.3 

20. 

-9.9 

14.4 

32.0 

42.4 

52.2 

52.2 

47.9 

49.6 

55.1 

60.6 

66.9 

73.7 

86.4 

77.6 

7.4 

25. 

-10.9 

12.9 

28.8 

41.2 

51.1 

55.3 

52.8 

47.9 

54.2 

57.3 

62.3 

69.3 

82.9 

75.8 

8.0 

30. 

-11.0 

10.9 

23.5 

35.0 

44.8 

50.9 

52.9 

47.4 

52.3 

53.3 

57.4 

64.2 

78.1 

75.8 

6.4 

35. 

-3.9 

9.0 

15.2 

22.9 

31.7 

42.5 

44.8 

45.7 

46.8 

48.9 

52.7 

59.0 

71.0 

75.4 

2.3 

40. 

-1.7 

3.6 

13.0 

14.9 

22.2 

31.2 

37.1 

40.0 

42.2 

43.0 

45.4 

52.3 

65.3 

73.6 

1.7 

45. 

1.0 

-1.3 

4.3 

3.5 

10.7 

18.6 

28.7 

32.7 

33.8 

34.5 

36.9 

44.0 

57.4 

70.4 

-3.7 

a 

0. 

5.4 

4.7 

6.5 

7.8 

9.3 

15.1 

LC/Q 
20.3 27.4 

31.4 

32.0 

28.7 

18.9 

8.6 

3.5 

4.3 

2. 

10.0 

9.6 

11.3 

14.6 

17.6 

23.5 

29.7 

37.2 

41.5 

40.3 

34.5 

20.9 

5.3 

3.6 

4.9 

4. 

13.6 

13.3 

17.9 

21.2 

26.3 

32.2 

38.6 

46.6 

51.8 

47.9 

39.0 

22.4 

3.2 

3.0 

4.0 

6. 

16.8 

16.2 

20.6 

25.6 

32.1 

38.7 

47.8 

55.0 

58.4 

51.8 

39.7 

23.9 

3.8 

2.8 

5.3 

8. 

19.7 

18.9 

21.6 

28.7 

36.1 

43.5 

54.1 

62.8 

64.4 

54.6 

38.9 

24.4 

4.5 

2.6 

5.6 

10. 

20.9 

21.1 

24.7 

30.8 

38.6 

48.1 

57.8 

64.8 

65.7 

55.0 

39.9 

26.3 

6.2 

2.0 

5.9 

12. 

24.2 

23.4 

27.2 

32.9 

41.8 

53.1 

61.3 

66.6 

61.9 

53.3 

41.9 

28.8 

8.0 

1.2 

6.2 

16. 

27.9 

29.0 

33.9 

40.7 

48.5 

61.4 

73.0 

67.8 

60.4 

51.3 

42.8 

31.8 

11.6 

1.1 

6.6 

20. 

35.5 

35.4 

39.8 

46.7 

54.0 

65.8 

73.3 

69.3 

62.1 

53.6 

45.4 

34.4 

15.2 

3.3 

7.1 

25. 

44.6 

43.2 

46.4 

52.8 

57.3 

63.6 

68.9 

71.1 

66.5 

58.7 

50.9 

39.5 

20.6 

7.6 

7.6 

30. 

53.4 

50.3 

51.8 

58.6 

61.5 

65.0 

65.1 

71.5 

68.8 

62.2 

55.4 

43.7 

25.7 

10.9 

8.2 

35. 

55.2 

56.9 

58.2 

59.5 

66.1 

68.0 

69.1 

70.9 

69.8 

63.3 

58.3 

45.7 

29.6 

14.4 

6.8 

40. 

61.4 

59.9 

62.8 

66.2 

67.5 

68.2 

69.9 

70.7 

67.3 

62.4 

59.1 

47.8 

34.2 

19.5 

9.0 

45. 

63.2 

64.0 

65.7 

66.7 

65.9 

68.0 

68.8 

67.9 

65.2 

60.5 

57.4 

47.7 

37.4 

23.1 

8.2 


0.0 

4.0 

8.0 

12.0 

16.0 

20.0 

25.0 

30.0 

35.0 

40.0 

45.0 

50.0 

60.0 

75.0 

90.0 
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Figure 2.2.- Wind tunnel data corrected for estimated support interference. 
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( b ) Moment components . 


Figure 2 . 2 .- Concluded . 
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Figure 2.4.- Corrected wind 














3. ANALYTICAL EXTRAPOLATION 

The symmetry of the MILVAN is analyzed in this section for equations to extrapolate the 
data in k 0 to the complete domain of interest for simulation: 

V = {(V>, a) : tp w E [-180°, 180°], a € [-90°, 90°]} 

It is useful to divide this domain into eight quadrants as shown and labeled in figure 3.1: 



* 

Figure 3.1.- Attitude domain for simulation models. 

An arbitrarily shaped load model can be turned on its side, upside down, and reversed in 
the wind tunnel to obtain data for all of D while only moving the model mount through angles 
in D 0 - For each orientation, the wind-axes force and moment components and the angles between 
wind and body axes are related to the components measured along the wind tunnel’s measuring 
axes and to the mount’s attitude angles. However, if the load is not arbitrarily shaped, but is 
a rectangular box with identical sides, then these various mountings are indistinguishable to the 
flow, and the force and moment components measured on the wind tunnel scales are necessarily 
identical in each orientation. In this case, the aerodynamics at any point in D can be defined 
from measurements at a corresponding point in D 0 . Finally, if a sufficient number of symmetry 
properties of the components are also assumed, then the aerodynamics at points within D 0 are 
related, and then the aerodynamics everywhere in D can be determined from those in a subregion 
of Do- It is found that the subregion with 0 < tp < 45° is such a region, but Z 0 is not. Nevertheless, 
the theory indicates relations which must be satisfied over D 0 by any simulation model derived 
from the measurements and approximation formulas, including the zero-lines of all components, 
and the values or behavior of each component along all boundaries of D 0 - 

We note that the ARC wind tunnel model did not have identical surfaces; the sides, bottom, 
and front were ribbed to represent strengthening ribs on the MILVAN while the top and rear were 
smooth. However, a limited amount of comparative data for an entirely smooth-surfaced box 
is given in reference 3 and indicates only small differences caused by the ribbing, and the same 
comparison is made in reference 6 with the same result. 


12 




3.1 Box Attitudes with Identical Wind Tunnel Measurements 
These attitudes can be obtained from the relations among four reference frames. First, 
the tunnel axes in which measurements are made are aligned with the tunnel longitudinal axis 
or free-stream direction, and the local vertical. Second, the mounting angles measured in the 
tunnel ( , a 0 ) define the mounting body axes, where the transformation from tunnel axes to the 
mounting axes is; 

Tb ot w o = •®'2( <a: o)^3(V , o) (3*1) 

The elementary rotation matrices, (Ei(o),i = 1,2,3) are defined in the nomenclature. The 
subcripts w 0 and b a indicate the tunnel and mounting axes, respectively, in the notation for 
vectors and transformations. 

Suppose the box is mounted with its body axes aligned with the mounting axes and oriented 
at some attitude (tp 0 , a 0 ) in D 0 , and the aerodynamic force and moment measured. If the box is 
now rotated about its longitudinal axis by any multiple of 90°, or about its vertical axis by any 
multiple of 180° , the box appears identical to the flow and the resulting measurements are all 
identical to those at (0 O , a 0 ). For all these orientations the corresponding transformation to body 
axes (indicated by the subscript 6) from mounting axes is 

Tb,b a = E3[m%)Ei(nir/2) m,n = 0,±1, ±2,... (3.2) 


The wind axes associated with the rotated box (indicated by the subscript, iui) are, by definition, 
oriented with the longitudinal axis in the free-stream direction and the vertical axis in the body 
vertical plane; consequently, the transformations from tunnel and body axes can be written, 
respectively, as 

T w>Wo = EMw) (3.3a) 

Tb,w = ^2(oti)£?3(V»i) (3.36) 


The angles <f> w ,oti and V’l are shown in the figure 3.2. Since i w and i Wo are parallel to Va, then 
tunnel and wind axes are related by a roll rotation about Va , which is denoted <f> w here. Since 
k w is in the body vertical plane then only two angles are required to locate wind axes relative to 
body axes, and these are taken as the usual angle of attack and (negative) sideslip angle as shown 
in the sketch. If the body axes are aligned with the mounting axes ((m,n)=(0,0)) then wind axes 
coincide with tunnel axes, <f> w = 0, and (V’ljO'i) = (V>o, <*o) 5 but this does not occur otherwise. 

Equations 3.1 to 3.3 can be combined to equate transformations between wind and body 
axes: 


Tb y w — Tb t b 0 Tb 0 ,w 0 T WotW 


from which the general equation for the attitudes (<f> w , V'ljOti) which have measurements identical 
to those at (V\>, oc 0 ) is; 


Ei(ai)Ez{^\)Ei{4> w ) = E 3 (mir)E 1 (nn/2)E 2 [a 0 )E 3 (xl> o ) m,n = 0,±1, ±2,... (3.4) 


Equation 3.4 yields nine scalar equations of which only five are independent. A convenient set 
of five equations used in this work is given in figure 3.3. For each case (m,n) two solutions for 
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WIND AND BODY AXES WIND AND TUNNEL AXES 


Figure 3.2.- Angles relating reference frames. 

(<f>u are obtained. Taking all cases together, only three yield distinct solutions and these 
are included in figure 3.3. For each point (fjf 0i a 0 ) in D 0 a single corresponding point in each 
quadrant is obtained, and, conversely, every point (V>i» «i) in V is related to a single point in D 0 
by equation 3.4. 

The force and moment components in wind axes are obtained from the tunnel measurements 
by the transformation T W)Wo . Using the notation D, Y, L, RM, PM and YM to denote the cus- 
tomary wind axes force and moment components, and recalling that drag and lift are positive 
along the negative longitudinal and vertical wind axes, then 

( ° \ 

— I Y cos <f>w L sm I 

.) yYsm^ + LcosKJ {l p ot0to) 

( RM \ 

PM cos <j> w + YM sin <f> w (3.56) 

-PM sin <(> w + YM cos <f> w J 



RM 

PM 

YM 



RM \ 

PM) =E 1 (<f> w ) 

YM I , 


3.2 Extrapolation from D 0 to V 

The various cases {(m, n)} in equations 3.4 and 3.5 suffice to relate the aerodynamics at any 
point in D to those at a corresponding point in D 0 . First, results for the case with 180° roll of 
the body axes relative to mounting angles ((m, ») = (0, 2)) relate the aerodynamics in D 2 and P3 
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EQUATIONS 


cos ai cos 0i = cos mit cos a 0 cos 0 O 

sin = — cosmir(sinn^s\na 0 costp 0 — cosn^sin0 o ) 

sin a\ cos 0i = sin n \ sin 0 O + cos nf sin a 0 cos 0 O 

cos <f> w cos 0i = cos m7r(sin n \ sin a 0 sin 0 O + cos nf cos 0 O ) 

sin <f> w cos 0i = cos mn sin n ^ cos a 0 


\ 


> m,n = 0, ±1, ±2,... 


Case: 

m 

0 

SOLUTIONS 

0 

1 


n 

2 

1 

0 

First 

Sol’n: 

0to 

0 

cot -1 (tan a 0 sin 0 O ) 

0 


01 

-180° + 0 O 

— sin -1 (sin a 0 cos 0 O ) 

-180° + 0, 


«1 

180° - a 0 

tan -1 (tan 0 O / cos a,,) 

~OL 0 

Second 

Sol’n: 

010 

180° 

— 180° + cot - 1 (tan a 0 sin 0 O ) 

180° 


01 

0 

1 

— 180° + sin - 1 (sin cx 0 cos 0 O ) 

0 

1 


<*1 

-«o 

—180° + tan -1 (tan 0 O / cos a 0 ) 

180° - a. 0 


Figure 3.3.- Box attitudes with identical measured aerodynamics. 
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with those in D\ and P 0 , respectively: 


D } 


( D \ 

Y 


-Y ' 

L 


-L 

RM 


RM 

PM 


-PM 

\YM) 

' 

(- 0 ,-a) 

\-ymJ 


( 0 , a) € Do, Di 


(3.6) 


That is, each component is either symmetric or antisymmetric about the origin of V along any 
radial line through the origin. 

Second, results for the case with 90° roll of the body axes ((m,n) = ( 0 , 1 )) define the aero- 
dynamics in V\ from those in D 0 \ 


where 


D \ 


D \ 

Y 


cos <p w Y — sin <p w L 

L 


sin<£„,Y + cos <p w L 

RM 


RM 

PM 


cos <p w PM + sin <p w YM 

\ym) 


— sin <p w PM + cos <p w Y M / 


(0,a) 6 D 0 


cot <p w = tan a sin ip = tan ai sin ip\ <p w £ [0, 90° ] 

sinV'i = — sin a cos ipi £ [— 90°, 0] 
tanaci = tant/>/cosa ai £ [0,90°] 


(3.7) 


Maps of (<p to, a i) over A> are given in figure 3.4. The wind axes roll angle varies from 90° 
on the boundaries a = 0 and ip = 0 to 0 on the boundary a = 90°. Referring to the figure 3.5, 
the points p 0 and pi are equivalent; the region R 0 in D 0 is equivalent to the region R\ in D\ 
(the left boundary of R\ is equivalent to the line a = 45° and is given by tan^i = — cosai); 
and their complement regions R 0 , are equivalent. Further, the lines a = 0 and ip = 0 are 
equivalent (this particular equivalence has been recognized frequently in the literature), and the 
line tp = 90° is equivalent to the single point (ipi,a i) = (0,90°) so that the aerodynamic 
components are necessarily fixed on this line independent of a. 

Third, the case with a reversal of box heading ((m, n) = ( 1 , 0 )) relates points in P 4 , D$, D& 
and D 7 to points in Pi , P 3 , D 2 and D 0 , respectively: 


D \ 


D \ 

Y 


Y 

L 


L 

RM 


RM 

PM 


PM 

\YMJ 

(—180 +<!>,-(*) 

\ym) 


\ip\ < 90° (3.8) 

) 
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0 15 30 45 60 75 90 0 15 30 45 60 75 90 


<P \l> 

(a) <t> w (4/,a): ROLL ANGLE FROM TUNNEL (b) EQUIVALENT POINTS IN^/D, 

TO WIND AXES 

Figure 3.4.- Equivalent points in Do and D\. 



Figure 3.5.- Equivalent regions of Do and D\. 

Equations 3.6 to 3.8 define the aerodynamics at any point in D from measurements at a 
corresponding point in D 0 - However, the available measurements are confined to the lower half 
of D 0 , and equation 3.4 does not relate pairs of points within D 0 nor provide a relation between 
every point in D and a point in R 0 . Additional symmetry properties of the box are considered 
next to obtain a relation between points in D 0 . 
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3.3 Symmetry Assumptions, Boundary Values, and the Equivalence of Measurements at 

Points within D 0 

Next, the following symmetry properties at the boundary between 
are assumed: 

Al. D is symmetric in ip about ip = 0 
A2. Y,YM are antisymmetric in ip about ip — Q 
A3. RM is antisymmetric in ip about ip = 0 

Properties Al and A2 are consistent with the available wind tunnel data and usually (but not al- 
ways) characterize flow around a body which is symmetric about its vertical plane. The symmetry 
of RM is less obvious a priori. It can be shown from equations 3.6 and 3.7 that a symmetric and 
an antisymmetric function have distinctly different properties at the boundaries of R 0 and within 
R 0 . The Ames roll-moment data are not consistent with either symmetry, but data from both 
reference 4 and 6 in appendices B and C show that it is antisymmetric and this is assumed here. 
It follows from equation 3.7 and property A2 above that L and PM are necessarily symmetric 
about this boundary, so that the complete symmetry properties about ip = 0 are: 


D 0 and Pi (at ip — 0) 


) 


(3.9) 


D \ 


D \ 

Y 


-Y 

L 


L 

RM 


-RM 

PM 


PM 

\YMJ 


\-ymJ 


(Vs <*) g Do 


(3.10) 


Symmetry about the D 0 — D 2 boundary (about a = 0) is obtained by combining equations 3.6 


and 3.10; 


D \ 


D \ 

Y 


Y 

L 


-L 

RM 


-RM 

PM 


-PM 

\ym) 

(0.-a) 

YM ) 


{ip, a) £ P 0 


(3.11) 


Naturally, the force and moment components are zero at the boundaries about which they are 
antisymmetric, assuming continuity of the static aerodynamics throughout D. 

The aerodynamics at pairs of points within 0 o can now be related. Referring to figure 3.6a, 
if {ip, a) is any point in Do, and (— ipi,<xi) is its equivalent point in V\ from equation 3.7, then a 
relation between the aerodynamics at {ip, a) and {ipi,ot\) is obtained by combining equation 3.7 
and 3.10; 


D \ 


D \ 

Y 


— cos <p w Y + sin^u,L 

L 


sin <p w Y + cos <p w L 

RM 


-RM 

PM 


cos <p w PM + sin (p w YM 

\ym) 


\ sin <p w PM — cos <p w YM / 


(VS <*) G Do 


(3.12) 
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where 


cot <f> w = tan a sin 0 = tanaisin0i 4> w G [0,90°] 
sinV’i = sin a cos 0 0i G [0,90°] 
tanai = tan 0/ cos a ai € [0,90°] 

Points related by equation 3.12 are termed “equivalent” hereafter. This relation is identical to 
that plotted in figure 3.5, except that values of 0i are now positive. Referring to figure 3.6b, 
the regions £2 and £2 are separated by the line of self-equivalent points (the line for which 
(0i, ai) = (0,a) in eq. 3.12) given by 

0* (a) = tan -1 (sin a) aG[0,90°] (3.13) 

Each point in £2 is equivalent to a point in Z 2 and conversely. Therefore, the aerodynamics at 
any point in D 0 (and any point in V) can be given from the aerodynamics at a corresponding point 
in % 2 ‘ Further, the left half of D 0 (the region, 0 < 45°) contains Z 2 and is more than sufficient 
to define the aerodynamics throughout D . 

The available measurements are in the lower half of V 0 (the region a < 45°). This region is 
equivalent to the region £1 (fig. 3.6c) whose interior boundary is equivalent to the line a = 45° 
and is given from equation 3.12 by 


0( 45 )(a)=tan l (cosa), a G [0,90°] (3-14) 

Unfortunately, this encloses only about one-third of the region a > 45° for which we have no 
measurements. 

The boundaries of R. 0 ,%i and Z 2 subdivide D 0 into six regions, and it is useful for the 
analysis to denote these separately as shown in figure 3.6d. The figure includes formal definitions 
of these regions for use in the calculations of the next section. The region of measurements, 
Ro, is subdivided into Si, $2 and S3. Each point in Si is equivalent to a point in S2 and the 
converse is true. Measurements at equivalent points within these two regions can be compared 
for their agreement with equation 3.12 and discrepancies are presumed to be due to measurement 
errors. Each point in S4 is equivalent to a point in S3 so that measurements in S3 will provide 
measurements of the aerodynamics at corresponding points in S4. Last, S5 and Sq are mutually 
equivalent, and the aerodynamics in these regions cannot be given from measurements in Z 0 
using equation 3.12. Thus, the symmetry properties represented by equations 3.4, 3.5, and 3.8 
are insufficient to define the aerodynamics everywhere in V from those in Z 0 . 

Although the aerodynamics in S5 and S& cannot be determined from measurements in % 0 , 
the theory and data provide values along the boundaries of S 5 and S& and an equivalence relation 
which must be satisfied by any approximate model of the aerodynamics for this region. 

Symmetry about the V 0 - boundary (about 0 = 90°) can be obtained by combining 
equation 3.8, which relates points in Dz and P 5 , with equation 3.6, which relates points in Vz and 


19 



(a) EQUIVALENT POINTS IN<D 0 



sin i//j = sina cos\p 
tan ay = tani///cosa 


(b) THE LINE OF SELF-EQUIVALENT POINTS 


\p* = tan - ^ (sina) a e [0, 90° ] 

q? 2 = {(^,a): [0, **(«)], a e [0,90°]} 


0 45 90 

* 

(c) <&,! REGION EQUIVALENT TO 9g 


^(45) - tan -1 (cosa) a e [0, 90° ] 

Qy = {<*.«> ^ e | 0 , ^ ( 45 , (a)],ae [ 0 , 90 °]} 


0 45 90 

'P 




(d) SUBDIVISION OF^D 0 FOR EXTRAPOLATION CALCULATIONS 

S-j = {{\p, a): \p e [0, \p* (a)] , a e [0,45°]} 

5 2 = {(^.a): \p e [\p* (a), ^< 45 >(a)] , a e [0,45°]} 

5 3 = {(*, a): * e [^ (45) (a), 90°] , a e [, 45°] } 

S 4 = (W.a): \p e [0, ^ 45| (a)],ae [45°, 90°]} 

S 5 = {(^,a): l\f/ (45, (a), f(«)],«e [45°. 90°]} 

S 6={(^,a): ^ e t^*(a), 90°] , a e [45°, 90°]} 



Figure 3.6.- Equivalent points and regions in Do- 
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(3.15) 
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(90— Aip , a) 
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(90+Aip,a) 


Y,L,PM and VM are antisymmetric in ^ about ip = 90° and, therefore, zero on this boundary, 
while D and RM are symmetric. As previously noted with equation 3.7, all components are 
constant on this boundary; since RM( 90°, 0°) = 0 (equation 3.11) we have 


Y = L = RM = PM = YM = 0 at ip = 90° (3.16a) 


D = D{ 90°, 0°) at ip = 90° (3.166) 

Last, symmetry about the boundary, a = 90°, can be obtained by combining equation 3.10, 
which relates points in V\ and P 0 , with the second solution of equation 3.4 for the case 
(m, n) = (1,0), which relates points in Pi to points in the quadrant (ip, a) € [0°, 90°] * [90°, 180°]; 
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(yp, 90— Aa) 
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(0,90+Aa) 


(3.17) 


The antisymmetric components, L,RM and PM are zero at a = 90°. An evaluation of equa- 
tion 3.12 on this boundary indicates that the symmetric components, D,Y and YM are either 
symmetric or antisymmetric in ip about ip = 45°: 


L = RM = PM = 0 



(45-A^O 


~ Y ] 

-Y **)<«+*» 


at a = 90° 


(3.18) 


Thus, Y and YM are zero at ip e {0,45°, 90°} along this boundary. 


3.4 Summary 

The geometric symmetry of the box results in a set of attitudes {(ipi,ai)} corresponding to 
rotations of the box in the tunnel for which the measured aerodynamic components are identical 
or differ only by a sign (eqs. 3.4, 3.5 and 3.9). This permits determination of the aerodynamics at 
any point in the set from measurements at one point, and at any point in P from measurements at 
a corresponding point in P 0 (eqs. 3.6, 3.7 and 3.8 or eqs. 3.8, 3.10 and 3.11), or at a corresponding 
point in the smaller region within P 0 given by 0 < ip < 45° (eqs. 3.8,3.10 and 3.12). 
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Measurements confined to the region within D 0 given by 0 < a < 45° do not suffice to define 
the aerodynamics throughout P; the aerodynamics at any point in R 0 or S4 (fig. 3.6) and its 
related region in each other quadrant can be obtained from such measurements, but not at any 
point in S5 or 56 (fig- 3.6) and their related regions. 

The exact aerodynamics and any model constructed from the available measurements must 
satisfy equation 3.12 and the zero boundary values specified by equations 3.10, 3.11, 3.16 and 
3.18; that is, in V 0 \ 


1. D satisfies equation 3.12. 1 

2. L and Y satisfy equation 3.12 and 
Y = 0 if 0 = 0 or 90°, 

L = 0 if a = 0 or 90°, or 0 = 90°. 

3. RM satisfies equation 3.12 and ► 

RM = 0 if 0 = 0 or 90°, or a = 0 or 90°, or tan 0 = sin a. 

4. PM and YM satisfy equation 3.12 and 
PM = 0 if a = 0 or 90° or 0 = 90°, and 
YM = 0 if 0 = 0 or 90°. 


(3.19) 


The measurements of reference 3 are reviewed in the next section for their degree of agreement 
with these conditions, revised as needed to satisfy equation 3.19, and then extrapolated consistent 
with equation 3.19 throughout V 0 to provide a comprehensive estimate of the box aerodynamics. 
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4. REVISION AND EXTRAPOLATION OF THE WIND TUNNEL DATA 

The object of this section is to estimate the MILVAN’s static aerodynamics throughout Do 
by fitting the Ames wind tunnel data in Ro (fig- 2.2) to the theoretical model of the previous 
section. This is done by: 1) revising the data tables in R 0 to satisfy the theoretical equivalence 
relations and zero force and moment lines (eqs. 3.12 and 3.19); 2) extending this fit to S4 using 
equation 3.12; and 3) defining approximation functions in Ss and S 6 which satisfy the known 
boundary values from equation 3.19 and the previous two steps, and the equivalence relations 
(eq. 3.12), and whose functional behavior is smooth at the boundary with S3 and S4 and simple 
within S5 and Se- The results are given as new data tables which span D 0 and fit the theoretical 
model. 

The Ames data in R 0 do not exactly satisfy equations 3.12 and 3.19. The discrepancies, which 
are large for some components, can result from measurement errors or differences between the 
assumptions underlying the theory and the test conditions. However, in all cases there is sufficient 
evidence from the Ames data or from independent wind tunnel tests (refs. 4-6) to show that the 
MILVAN’s aerodynamics agree closely with the theory. Therefore, all significant discrepancies 
are assumed to be due to measurement errors and the estimate in R 0 is obtained by revising the 
data to fit the theoretical model. These revisions are small for drag, sideforce, and yaw moment, 
moderate for lift, and the roll- and pitching- moment data are entirely discarded and replaced 
by estimates obtained from independent data (roll) or from the remaining data and the theory 
(pitch). The resulting estimates are all compared with the available independent wind tunnel 
data, and good to fair agreement is obtained. 

An analysis sufficient to identify the systematic and random tunnel errors in measuring forces, 
moments, and attitude angles or the data processing errors, which might account for the principal 
differences between the Ames data and the theory, is beyond our scope. The original tunnel records 
are no longer available for such an analysis nor was there any opportunity for further tunnel tests. 
In several cases, systematic measurement errors consistent with the observed discrepancies along 
the theoretical zero lines are noted as an aid to selecting plausible global behavior of the correction 
functions, but very little can be concluded from the data alone regarding actual error sources. 

4.1 Drag 

The drag data can be tested for their agreement with equation 3.12 over the regions $1 and 
52 by calculating the difference between measurements at equivalent points: 

6D(if),a) = DC(ip,a) - DC(i() 1 ,on) (ip,a) € Si or $2 (4.1) 

where: sin tpi = sin a cos rjs 

tan ai = tan if)/ cos rf) 

and DC is the corrected wind tunnel data from figure 2.2. These differences (fig. 4.1) are small 
everywhere in 5i and 52 with a maximum value corresponding to 5% of the measured drag. 

An estimate of the drag in 5i and 5 2 which satisfies equation 3.12 is obtained by averaging 
the data at equivalent points: 

D(tjj, a) = 0.5 (DC(0, a) + £>C(0i, a 1 )) (V>, a) G 5i or S 2 (4.2) 
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-A. 

This revises the drag data only slightly (DC — D = 0.561?) . 

To examine the appearance of measurement errors in these revisions first express the data, 
DC, as the sum of the actual drag, D, and the drag measurement error, DC, and then obtain: 


D(tp,a) = D(rJ;,a) +0.5(DC(iP,a) + DCtyuoa)) ' 
6D(1>, a) = DC(i>, a) - DC(ih,ct x ) 


(4.3) 


Thus, D retains the average measurement error at the two equivalent points while 6D is their 
difference. This difference is also a lower bound on measurement errors; 


0.5\6D(x/j,a)\ < max{\DC (0,a)|, \DC(ipx,oti)\} (4.4) 


Thus, the measurement errors are at least as large as 6D/2 , and small values of 6 D admit 
small errors, but do not preclude a large table bias or average error. However, a comparison of 
the drag estimate, D, with independent data from references 4-6 in figure 4.2 show only small 
differences which indicate the probable size of the measurement errors and the effects of differences 
in surface details of the wind tunnel models. A more comprehensive comparison with the data 
from reference 6 is made in appendix B (fig. B3) and shows similar good agreement throughout £o* 
Further, the data from references 4-6 were taken at negative heading so that these comparisons 
also confirm the symmetry of drag with ip about ip = 0. 

An estimate of the drag in 53 and 54 which satisfies equation 3.12 is obtained by averaging 
the drag measurements at ip = 90°, denoted Dgo, and then 

(DC(iP,a) (xl>, a) € 5 3 ip ^90° 

Dty, a) = < D go ip = 90° (4.5) 

\DC(ip 1 , 01 ) (ip,a)eS 4 

Note that the drag measurements (fly. 2.2) are very nearly constant at ip = 90°, as expected from 
equation 3.12, and averaging to impose this boundary condition revises the data negligibly. 

Last, drag must be extrapolated over S 5 and Sq subject to equation 3.12. Consider, first, an 
extrapolation based on the theory that the drag of bluff bodies varies principally in proportion to 
the cross-sectional area of separated flow. This area is estimated here as the box area projected 
on the plane perpendicular to the air velocity vector: 


Ap/S = O.4|cos0cosa| + IcosV'sinal + ) sin 0] (4.6) 

where S is the side area, 160/f 2 . A plot of the measurements versus A p (fig. 4.3) shows good 
correlation at low values of A p but much larger spread in the data at high values, and is nonlinear. 
However, when plotted versus the angle between the velocity vector and the box longitudinal axis; 

£ = (i Wo ,i bl ) = cos - 1 (cos cos a) (4.7) 
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better correlation is obtained (fig. 4.3b). Thus, the drag measurements can be fitted with a 
monotonic function versus £ with moderate discrepancies, and this fit can be extrapolated to a 
global model for drag in only one variable; 

D{rp, a) = D(£(V», a)) (ip, a) e Do 

where D(£) is the fitted function. It can be shown that both A p and £ are invariant at equivalent 
points (A p (ip,a) = A p (ip i,<*i) and £(ip, a) = £(^i,ai)) so that global models obtained by fitting 
the data to either variable would satisfy equation 3.12. 

The following estimate of drag over Ss and $e is adopted based on the measured drag versus 
£ along the boundary a = 45° ; 

D(ip,a) = Z?(V>(^,a),45) (tp,a) € S$ or Sq (4.8) 


where 

ip(ip,a) = cos -1 (v^ cos £(Vb<*)) 

Here, D(ip, 45), £ are defined in equations 4.5 and 4.7 above, and rp is an inversion of equation 4.7 
at a = 45°. This extrapolation satisfies equation 3.12 since it is a function of £ only. Further, it 

A A 

is continuous and “smooth” at the boundaries with S 3 and S 4 (that is, D in eq. 4.8 matches D in 
eq. 4.5 at these boundaries, and the functional behavior of D does not change abruptly at these 
boundaries as can be seen in the contour plot in fig. 4.4), and the extrapolation assigns values 
which are within the range of the measured drag values. 

The drag estimate defined by equations 4.2, 4.5 and 4.8 is given in figures 4.4 and 4.5. This 
estimate matches the measured behavior as closely as possible while satisfying equation 3.12. 

4.2 Side force and lift 

The side-force data and plots of figures 2.2 and 2.3 show small discrepancies from the ex- 
pected zero crossings at the the boundaries, rp = 0,90°. These discrepancies are assumed to be 
measurement errors which are a sum of random and systematic errors in measuring angles and 
force. Plots of the data for these boundaries indicate they are systematic and qualitatively con- 
sistent with a bias of 1° — 2° in measuring yaw angle, Arp. This results in a side- force bias at 
the end points of the form; b y = dY/dipAip. We assume these biases are continuous and affect 
measurements at other values of ip, so a correction function, given as a linear combination of the 
end point discrepancies, is defined; 

M*. <*) = (!- 55) rc ( 0 , a) + £ VC( 90, «) (4.9) 

and then 

YC 1 = YC-b y 

This revision nulls Y C at rp = 0° and 90° and otherwise varies linearly in heading between its 
values at these end points. This variation approximately follows the variation of BY j dip. A 
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comparison of the revised measurements with independent data from references 4-6 in figure 4.2 
shows good numerical agreement at a = 0, and a more comprehensive comparison in appendix B 
(fig. B4) shows good agreement throughout Rq. This implies that the data from reference 3 
do not contain a large systematic error at intermediate heading angles that was missed in the 
end-point data used to define the revision in equation 4.9. 

The lift data in figures 2.2 and 2.3 show good agreement with the expected zero value at 
xp = 90° , but agreement with the expected zero value along the a = 0 boundary is poor; there the 
measured lift reaches 32 Q lb at intermediate values of heading angle. This variation of measured 
lift with heading at a = 0 is seen to be systematic and present consistently at a up to 30° in 
figure 2.3. It was noted that the box surfaces are not identical as assumed in the extrapolation 
theory of the previous section; the bottom, side, and front of the wind tunnel model were cor- 
rugated while the top and rear were smooth. This might result in a zero lift line different from 
a — 0. However, reference 3 contains some data for a smooth-surfaced box which was compared 
with the data of figure 2.2 and shows the same variation of lift with xp at a = 0 and generally 
small differences from the corrugated box lift over the region for which data are given. That is, 
the effect of the surface corrugation on lift is much smaller than the maximum lift measured at 
a = 0. Further, data from reference 6 is available for both smooth and corrugated containers 
and these also show only small differences in lift, but good agreement with the expected zero-lift 
line at a = 0 (appendix B, fig. Bl). Finally, it is noted that the questionable lift variation in 
the Ames data was consistent with a sizeable systematic error in measuring angle of attack, Aa, 
that is; the variation of lift at a = 0 approximately follows the variation of dL/da with yaw. A 
correction was defined from the nonzero boundary values as: 


M0,«) = 95 LC ( 90, a) + |iC(0,O) - i tC(90, 0)|(1 - £) 
LC1 — LC — b L 


(4.10) 


This revision nulls LC at the boundaries xp = 90° and a = 0. The second term varies linearly in a 
and approximately follows the variation of the lift slope, which declines monotonically with a to 
small values at a = 45°. The resulting table of values of LC1/Q , contained small negative values, 
under 5 ft 2 , in the region xp > 60°; these were assumed to be due to errors in estimating the 
systematic errors rather than present in the actual lift function, and were faired out of LC 1 for 
the following calculations. The discrepancies from theoretical boundary values noted previously 
and the independent wind tunnel data suggest that lift is less accurately measured than side force 
in the Ames data. 

The discrepancies between the revised data and values calculated from the data at the equiv- 
alent point are: 

(YC1\ _(-cos<p w sin^\/yCl\ , , 

\ LC 1 J^a) COB A /(*,«.) { } 

where (xp, a) € Si, S 2 ; Y Cl, LC 1 are defined in equations 4.9, 4.10; and (xpi, <* 1 , <p w ) are functions 
of (xp, a) defined here and everywhere in this section by equation 3.12. A table of lift discrep- 
ancies is included in figure 4.1. The results show moderate discrepancies under 15Q at nearly 
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all gridpoints in Si and S 2 - Plots of the two lift functions compared in equation 4.11 are given 
in figure 4.1c and show adequate agreement in curve shapes. The agreement of these data with 
equation 3.12 is somewhat poorer than for drag. 

To abbreviate the notation hereafter, let /cl,/, etc., be the two-dimensional force vectors 
(YC 1 ,LC 1 ) T , ( Y,L ) T , etc., and let U be the coefficient matrix for these components in equa- 
tion 3.12 (repeated in eq. 4.11). 

An estimate, /, of these components which satisfies equation 3.19 can be given from /cl by; 
1 ) averaging /cl at equivalent points in Si and S 2 to impose equation 3.12, 2) using equation 3.12 
to define / in S 4 from /cl in S 3 , and by 3) interpolation of / in S 5 , Sq, between its known values 
along the boundaries of this region, assuming linearity in a, followed by averaging at equivalent 
points in S 5 and Sq to impose equation 3.12; 


/(V>, a) 


' 0.5(fcl(ip,a) + U fc l(V»i,ai)) 

. fcl(rp,a) 

‘ Ufeltyuai) 
lo.5(//(^,a) + tr//(tfi,ai)) 


(V>,a) € Si, $2 
(0,«) € S 3 
(ip, a) € $4 
(ip, a) e Ss, $6 


where fl is the linear function defined in S 5 and Sq by; 




(ip,a) e S 8 , Sq 


ao 


f 45 

\ cos - 1 (tan(^)) 


ip > ip* (4 5) 
ip < ip*( 45) 


(4.12) 


Here, ip* is defined by equation 3.13, ip*( 45°) corresponds to the meeting point of all regions 
Si — Se seen in figure 3.3, and {ao(ip); 0 < ip < 90°} is the lower boundary of S 5 and Sq. 

The estimate in Si and S 2 is obtained by averaging /cl with the value calculated from /cl 
at the equivalent point. The appearance of measurement errors in this estimate can be obtained 
by representing the revised data, /cl, as the sum of the exact force, /, and the error, /cl, and 
then, assuming equation 3.12 is exact, / can be written as 

f(ip,ci) = f{ip,ci) + Q.5(fcl{ip,a) + Ufcl{ipi,ai)) (ip, a) € Si, S 2 (4.13) 

where / contains the average of the errors at the two points. Note that U is orthogonal and does 
not change the magnitude of fcl(ip\, <* 1 ), but can distribute this magnitude in any ratio to the 
components Y and L. 

Alternatively, the estimate in Si and S 2 can be given solely from the side force data by solving 
equation 3.12 for L(ip,a) in terms of Y(ip,a), and y^i^i); 

Y\iP,a) = YC 1(0, a) 

L'(ip, a) = cot <p w Y Cl(ip, a) + esc <p w Y Cl(ip\ , ai) 
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however, this lift estimate is ill-conditioned at small <f> w , and the accuracy of V is degraded from 
the Y-measurement accuracy. An estimate based solely on the lift measurements can also be given 
analogously with a similar effect on estimation accuracy. 

A comparison of the side force estimate at a = 0 with independent data is included in 
figure 4.2, and a more comprehensive comparison of / with the data from reference 6 over the 
region £ 0 is given in appendix B (fig. B4). Good agreement is obtained for sideforce (differences 
are below 10Q at nearly all grid points). Lift differences are larger, in the range of 15Q to 30Q 
over much of £ 0? but these differences are the same size as the disagreement of the lift and side 
force data from reference 6 with equation 3.12, and are probably due primarily to lift measurement 
errors in the data from reference 6. 

-a, y-v 

The extrapolation of / into $ 5 and Sq is carried out as follows. First, / is assumed zero on 
the upper boundary of S5 and 56 (at a = 90°). For lift, this value is required by equation 3.18. 
For side force, this value is consistent with equation 3.18, which requires Y = 0 at the three points 
1 p € {0,45°, 90°}, and the expectation that Y is small, if not zero, on this boundary. Second, 
/ is known throughout Ro and S4 from the measurements; this data can be plotted versus a 
and the resulting curve extrapolated to the known boundary value at a = 90° such that the 
extrapolation is both continuous and smooth at the lower boundary of S5 and Sq. This first 
extrapolation does not necessarily satisfy equation 3.12, but equation 3.12 can be imposed by 
averaging the function with the value calculated from the function at the equivalent point. In the 
present computations / is extrapolated linearly in a and this was found sufficiently smooth and 
computationally convenient. 

ys. A. 

The estimates, Y and L, defined by equations 4.9, 4.10 and 4.12 are given in figures 4.4 
and 4.5. The maximum lift of the MILVAN occurs outside the domain of the original measurements 
and is predicted from the side-force data. Maximum values exceed 80Q lb at a = 60° and low 
values of ip. The contour plots (fig. 4.4) show smoothly varying functional behavior throughout 
Do, and a comparison with the measurements plotted in figure 2.4 indicates that little qualitative 
change has been made to functional behavior over Zo, while numerical changes are small for side 
force and more significant for lift, principally near a = 0. 


4.4 Roll moment 

The roll moment data from references 4-6 (appendices B and C) indicate that RM is an 
antisymmetric function and, in £ 0 , it satisfies; 

RM{ip,a) = -RM{tp u oti) 

RM = 0 if ip = 0, or ip = 90°, or a = 0 

However, the Ames roll moment data (figs. 2. 2-2.4) have large discrepancies from these properties. 
Large magnitudes occur along all boundaries where zero is expected, and the systematic trend 
versus ip to large negative values which occurs at a = 0 is present at all a. 

To examine the symmetry of the Ames data further, the table can be separated over Si and 
S 2 into components which are symmetric and antisymmetric with respect to equation 3.12; 

RMCS{xP,<x) =0.5{RMC{tp,a) + RMC{rp uai )) 

RMCA(ip,a) = 0.5 {RMC{ip,a) - RAfC(^i,a x )) 


} 


(4.14)) 
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where RMCS and RMC A are the symmetric and antisymmetric components of the data, RMC. 
Any arbitrary table of numbers can be separated in this way. 

Assuming the exact function, RM , is antisymmetric and satisfies equation 3.12, the separation 
can be related to RM and the measurement errors, RMC , as follows 

RMCS {ip, a) = 0.5 {RMC {ip, a) + RMC{ip x ,a x )) 

_ _ ► (4.15) 

RMCA{ip, a) = RM{ip, a) + 0.5 {RMC{ip, a) - RMC{iP x , e*i)) , 

Here, RMCS contains the average error, and RMC A contains RM and the antisymmetric com- 
ponent of the measurement errors. Note that a table bias {RMC{ip,a) = RMC{ip\,a x )) or 
nearly symmetric errors are self-canceling in RMC A. 

The results (fig. 4.6) show that the measured antisymmetric function, RMC A, has large 
nonzero values to 100Q ft-lb along the boundaries, ip = 0 and a = 0, where zero is expected. 
Since these are among the largest values in the RMC A table then RMC A is predominantly the 
antisymmetric component of the measurement errors and any estimate of RM from RMC A would 
be completely uncertain, except for the zero values imposed by the theory. The average errors, 
RMCS, are also large, reaching 70 Q ft-lb, and exceed the measured antisymmetric function at 
many points in the table. The average error is a lower bound on the measurement errors; 

\RMCS{ip,a)\ < max{\RMC{ip,a)\, \RMC{ipx,<xi)\) (4.16) 

Thus, the roll-moment data from reference 3 appear to be dominated by measurement errors. 
Finally, we adopt the estimate 

RM{ip,a) = 0 {ip, a) €. Do (4.17) 

since 1) the data from reference 3 are uncertain, and 2) the data from reference 6 (appendix B) 
show much better agreement with the theory and indicate that roll moment is small throughout 
£ 0 ; therefore it will have little or no effect on load equilibrium and motion compared to the roll 
restoring torque exerted by a suspension. 

4.5 Yaw and pitching moments 

The yaw-moment data shows good agreement with its expected antisymmetry in ip about 
ip = 0 and 90° (fig. 2.3). Numerical values (fig. 2.2) show moderate discrepancies from their 
expected zero values at these yaw angles. These appear to be systematic and were corrected 
similarly to the side-force data; the measurement errors were estimated as 

bruit, a) = (1 - ^)YMC{0,a) + ^YMC{ 90, a) (t/>,a) € So (4.18) 

and then 

YMCl = YMC - b YM 
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This correction nulls Y MC a.t xp = 0 and 90° and otherwise varies linearly between the apparent 
end-point errors. A comparison with data from references 4-6 at a = 0 is included in figure 4.2 
and shows good agreement in functional variation, with moderate numerical discrepancies (be- 
low lOOQ ft- lb) among the three sources. A more comprehensive comparison with data from 
reference 6 in appendix B, figure B7, shows a similar degree of agreement throughout Ro- 

Pitching moment is predicted by the theory of section 3 to be antisymmetric in a about 
a = 0, and in xp about xp = 90°. The data in figures 2.2 and 2.3 show systematic departures from 
the predicted antisymmetries and large discrepancies from the predicted zero values along these 
boundaries. These discrepancies are also the largest measured values in the data table, and they 
are not due to the surface corrugations since the data for smooth and corrugated-surface boxes 
in reference 3 differ by amounts which are much smaller than these discrepancies. Independent 
data from references 4-6 in appendices B and C show small to moderate discrepancies from the 
expected zero-PM lines, and moderate discrepancies from the predicted relation with the YM 
data when tested throughout the region Ro- 

Consequently it is assumed that the discrepancies of the PM data in reference 3 from the 
theoretically predicted zero values are due to measurement errors and that these dominate the 
data at heading above 25°. Therefore, the PM data are discarded, and we estimate PM and 
YM entirely from the YM data and its extrapolation to the region a > 45° as follows: 

YM(xli,a) | FA/Cl(t/>,45)(90 - a)/45 (iM) S S 4 


PM(xp,a ) = cot xp w YM(xp,a) + csc<p w YM(xpi,a\) 


{rp,a) e Ro, S 4 


(4.19) 



/ YMCl \ / cos <p w 

\PMCl) { ^ a) + \ sm<p w 


sin <p w \ ( YMCl \ 

-cos <p w ) \PMCl) ai) 


where, in S 5 and Sq 

YMCl(xp,a) = YAfCl(V>,45)(90 - a)/45 
PMCl(xp,a) = g(xp,a ) = faired extrapolation function 


(ip, a) 6 S 5 ,S 6 


The estimate of PM in Ro and S 4 is obtained from YMCl at equivalent points using a 
relation obtained by rearranging equation 3.12. This relation is singular at <p w = 0, and otherwise 
can magnify the YM-estimation errors by the factor esc (p w so that PM is necessarily less accurate 
than YM. However, the estimate given by equation 4.19 for the regions Ro and S 4 is satisfactory 
since (p w < 45° there and the error magnification is under 1.5. The resulting pitching-moment 
estimate, PM, is about as accurate as the yaw- moment estimate. This estimate differs completely 
from the original data in functional behavior (compare PM in figs. 2.4 and 4.4), but a comparison 
with independent pitching-moment data from reference 6 over the region R 0 (appendix B, fig. B7) 
shows the same degree of agreement as seen in the comparison of yaw moment estimates in 
figures 4.2 and B7. 
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The estimate in $5 and Se uses the same extrapolation scheme employed for lift and side force 
in this region. Along the upper boundary, (a = 90°), PM is zero (eq. 3.18) and YM is assumed 
zero based on its known zero value at three points on this boundary (eq. 3.18) and the expectation 
that its value is small elsewhere on this boundary, even if not zero. In addition, YM and PM 
are known throughout £0 and S 4 from equation 4.19(a), and can be plotted versus a and then 
extrapolated across S 5 and ,Se to the known boundary values at a = 90°. This extrapolation, 
denoted YMC 1 and PMC 1 in equation 4.19(b), is made with little knowledge of the actual 
functional behavior, such as local extremes, in S 5 and S 6 . Therefore, a curve is faired between 
the known boundary values subject to; a) continuity and smoothness (continuous derivatives) 
at the lower boundary, and b) the minimum number of curve reversals needed to connect the 
boundary values. In the application of this conceptual scheme, Y M was extrapolated linearly in 
a while PM, which rises toward a maximum near a = 45°, required some nonlinearity in the 
faired extrapolation function (denoted as g(0, a) in equation 4.19 without further definition and 
implemented as a tabulated function obtained from manually faired graphs in the actual work) to 
obtain smoothness at the boundary. The resulting initial extrapolation, YMCl and PMCl, was 
then averaged at equivalent points in S 5 and Se to impose equation 3.12. The results are included 
in figures 4.4 and 4.5. 
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(c) Lift comparison. 

Figure 4.1.- Differences between data at equivalent points; force components. 
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( a ) Force components . 

Figure 4 . 5 .- Estimated MILVAN aerodynamics . 
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Figure 4.6.- Roll moment data separated into symmetric and antisymmetric components. 
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5. DISCUSSION AND CONCLUSIONS 
Thus far, the results consist of a table of wind-axes force and moment components in the 
base domain, Do, (fig. 4.5) together with formulas in section 3 (eqs. 3.8, 3.10 and 3.11) for 
extrapolation to the complete domain of interest for simulations. These formulas are collected 
together in figure 5.1. For any point (ip, a) the force and moment components are given from 
the data at a corresponding point (ipot&o) in A> by a relation which depends on the quadrant of 
(V>, a). The form of this relation together with the correspondence of points and the coefficient 
matrix for each quadrant are listed in figure 5.1. The coefficent matrices are all diagonal and 
represent rotations of 180° about one or another axis. A simulation routine is readily written 
from this information. 

A body axis formulation is often preferred in simulations. The required data tables are not 
given here but can be generated routinely from the data in figure 4.5 by the transformation 




(5.1) 


where 


( cos a cos V > 

-sint/) 

sin a cos ip 


cos a sin ip 

cosip 

sin a sin ip 


— sin a A 



cos a J 


and / refers to either the force or moment vector. The corresponding extrapolation formulas can 
be obtained by applying equation 5.1 to the extrapolation formulas for wind axes. The results 
are included in figure 5.1; these have the same form as for wind axes, but the coefficient matrices 
are modified for quadrants with \ip\> 90°. 

In addition, sideslip angle is customarily used in simulations in place of heading relative to 
wind axes. This can be accommodated by revising signs in the data tables for either wind or body 
axes in accordance with equation 3.10 and then equations 3.8, 3.10 and 3.11 and the formulas of 
figure 5.1, with /? replacing ip, govern the extrapolation to all quadrants. Another possibility is to 
formulate the model in vector form without special notation for each scalar component; for this, 
the signs of D and L in the data table are changed, and then the formulas of figure 5.1 apply 
unchanged for both wind and body axes; e.g., for wind axes 


( FA W \ [FA W \ 

\ MA J (M \ ma J ( ^> 

where FA and MA are the aerodynamic force and moment vectors. 

Various indications of the accuracy of the data and the present estimate are noted in the text 
and summarized next. First, differences in repeated measurements indicate the presence of random 
measurement errors and unsteady flow effects not averaged out in the data processing. Second, 
discrepancies of the data from the zero lines and from the relations between the data at equivalent 
points predicted by the theory of section 3 reveal the presence of systematic errors because of 
sensors, test conditions, and data processing, and from differences between the assumptions of 
the theory and the test conditions. Third, differences in data and results from the independent 
data of reference 6 reflect errors in the data from both sources, and provide an independent 
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indication of the systematic errors. The data from reference 6 do not show the large-to-gross 
discrepancies from theory found in the Ames data for L, RM and PM, but are not otherwise in 
closer agreement with theory. The effects of nonidentical model surfaces are present in the second 
and third comparisons, but these effects were measured in both wind tunnel studies and were 
found to be small. 

The drag data show good repeatability and good agreement with the expected symmetry at 
equivalent points in the region for which this can be tested (in Si and S 2 ), and the drag estimate 
agrees closely with the data from reference 6. The side-force data show good repeatability and 
agreement with the expected zero lines and with the data from reference 6. The lift data show 
good repeatability, good agreement with the zero-line at ip = 90° but significant discrepancy from 
the zero line at a = 0° for intermediate heading and, when these disagreements are corrected, the 
results show good agreement with the theoretical relation of lift and sideforce measurements at 
equivalent points. The lift estimate shows only fair agreement with the lift data from reference 6, 
but the difference is likely dominated by inaccuracies in the lift data from reference 6. The roll 
data show poor repeatability and no agreement with any symmetry property or with the data 
from reference 6, and the pitching moment data show poor repeatability and gross disagreement 
from the expected zero lines at a = 0, and ip = 90°. In both cases, the largest values in the data 
table are presumed measurement errors and these data are considered unusable. Finally, the yaw 
moment data show good repeatability and agree well with the expected zero lines at ip = 0 or 90°, 
and the estimates of both PM and Y M obtained from the Y M data show fair agreement with 
the data from reference 6. In summary, it appears from these comparisons that the accuracy of 
the estimate is good for drag and side force and fair for the remaining components in the regions 
for which the comparisons can be made. 

A more definitive assessment of the inaccuracies in the data from reference 3 and a more 
accurate correction of that data than the revisions of section 4 requires a study and identification 
of the principal underlying tunnel and data processing errors. Unfortunately, the possibility of 
such a study is limited since the measurement and data processing records have been discarded. 
A review of the equations relating the six wind axes components of the aerodynamics to the six 
tunnel scale outputs for the Ames 7 x 10 ft tunnel indicated that the gross trends of RM and 
PM versus ip seen in figure 2.3 (which are also the systematic discrepancies of these data from 
the theory) Eire consistent with a significant error in the vertical coordinate of the model support 
point in the data processing. This results in RM and PM errors proportional to side force and 
drag, respectively, and has no effect on any other component. A correction of the RM and PM 
data based on this model of one error source reduced the discrepancies, but large discrepancies 
remained apparently because of significant additional systematic and random errors. 

Factors which limit the accuracy of the present estimate of the MILVAN’s static aerodynamics 
are; 1) random and unidentified systematic tunnel errors in the available data, 2) the sparseness of 
the heading-angle measurement grid in the region 45° < ip < 90° , which may result in significant 
interpolation error in generating Y and Y M from the data table, as well as a lack of independent 
comparison data to confirm the estimate in this region, and 3) the lack of data for the region 
« > 45°. Model uncertainty is greatest in the region a > 45°. The analysis of symmetry shows 
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that the aerodynamics in most of this region are unrelated to the aerodynamics in the region of 
available measurements. The present model is defined in the subregion S4 from the measurements 
in S 3 , but this could only be done for the force components, and otherwise, the model for a > 45° is 
defined by fairing curves between theoretical or measured boundary values. Thus, it appears that 
the present global model is subject to small to moderate uncertainties depending on component 
and region of the domain, and that improvements in accuracy require new measurements with 
improvements in measurement accuracy, tunnel error corrections, and measurement grid density 
and domain. 

Finally, the following conclusions are drawn: 

1 . The symmetry of a rectangular box permits determination of its static aerodynamics at any 
attitude from its aerodynamics at a corresponding point in a restricted attitude domain in Vo- 
The smallest such domain is 

{(^,a); 0 < ip < tan -1 (sina:), 0 < a < 90°} 

This is included in the region 0 < t(> < 45° of Vo, which is more than sufficient to define the 
aerodynamics globally. However, the measurement grid of the Ames data encloses the region 
0 < a < 45° of V 0 . This suffices to define the aerodynamics everywhere except in 

{(0,a); tan -1 (cos a) < if> < 90°, 45° < a < 90°} 

and its equivalent regions in other quadrants of V\ that is, except in about two-thirds of the region 
of V with |a| > 45°. 

2. For all components, zero lines and/or relationships between data at equivalent points are 
determined by the symmetry of the box. This feature permits comprehensive testing of the wind 
tunnel data for consistency with these properties throughout the attitude domain. Discrepancies 
indicate the presence and size of measurement errors in the data. 

3. For the Ames data, discrepancies are small for drag, side force, and yaw moment, large for 
lift, and gross for roll and pitching moment. Limited data from a University of Maryland study 
agree more consistently with the theory, and discrepancies range from small to moderate in size, 
depending on the component. In addition, repeated measurements in the the Ames data show 
small differences for the force components and yaw moment, and large differences for the RM 
and PM data. 

4. The accuracy of the global model estimated from the Ames data varies with component and 
over the domain. Accuracy is judged to decrease by component in the order drag, side-force, and 
yaw moment, lift, pitching moment, and roll moment. Accuracy over the domain is best at small 
angles, |a|, \xp\ < 45° where most of the measurements are concentrated, and poorer elsewhere 
because of larger yaw measurement grid intervals (10° to 15°), or the absence of measurements. 

5. The existing data provide adequate accuracy for the most significant components of the aero- 
dynamics in the region surrounding the minimum drag attitude of the MILVAN and suffice for a 
global model of moderate accuracy. Gains in the accuracy over that obtained here require new 
measurements with improved measurement grids. 
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S t = diag{l, —1, 1, —1, 1, — 1} = diag{-E 2 {x),E 2 {x)} 

= dtag{l, 1, -1, -1, -1, 1} = diag{-E 3 (x),E 9 (x)} 

5 3 = diag{ 1, -1, -1, 1, -1, -1} = diag{Ei(ir), Ei(w)} 

5 4 = dtag{-l, 1, 1, 1, -1, -1} = d»aff{-JSi(ir), £i(w)} 

5 5 = diag{-l, 1, -1, -1, 1, -1} = dtasf{£*(ir), £? 2 (*)} 

5e = dtag{- 1, -1, -1, 1, 1, 1} = diag{-1, 1} 

S 7 = diag{- 1, -1, 1, -1, -1, 1} = dtag{E 3 (ir), £ 3 (ir)} 

Figure 5.1.- Equations for the extrapolation of data in P 0 to any point in P. 
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APPENDIX A.-AMES RESEARCH CENTER WIND TUNNEL DATA 

The reference table of wind tunnel data used in the text is obtained by combining the data 
provided in reference 3 for all points in the reference domain, 

Zo = {U>,a)i V'€[0°,90°], a 6 [0°, 90°]} (Al) 

This combining can be done in two steps as described below. 

Al. Basic Measurements 

First, reference 3 provides two tables of measurements of the force and moment components 
in wind axes, and these must be combined. The two tables are obtained from measurement 
sequences made with fixed angle of attack (yaw-sweep data) and fixed yaw angle (pitch-sweep 
data), respectively. These data are quoted in figure Al for all points which fall in k 0 . As seen 
from the grids in figure Al, the two tables together cover all of k 0 and have a common region 
over which two sets of measurements are available, but at nonidentical grid points. 

These data can be combined by; 1) defining a reference grid as the union of the two grids, 
and 2) averaging the two tables in the common region, using table interpolation if a reference 
grid point is absent from one or both tables. This treatment uses all available data, halves the 
discrepancy between the two tables, and the combined grid provides better function definition 
than either individual grid. 

Errors in the data are assumed to be a superposition of repeatable systematic errors which can 
vary with attitude, and random errors which are independent from measurement to measurement 
and are assumed to have similar statistics for both tables; that is, the error in a measurement of 
/ at (ip, a) can be written as: 

f(rl>,a) = f 0 (t!>,a) + r) (A2) 

where /, /„, r? are the total, systematic, and random errors in the data, respectively. These errors 
can arise from multiple causes originating in the sensors, tunnel flow, test apparatus, and data 
processing, but the principal sources cannot be identified here since the test conditions are no 
longer accessible and the analysis is confined to the measurements. In the tunnel tests, the scale 
readings varied continually because of unsteady flow effects, and the measured static aerodynamics 
were obtained by averaging numerous readings. Similar unsteady effects and their treatment by 
averaging occurred in all investigations; the magnitudes of these effects have not been reported, 
but we assume they are small, and that the static aerodynamics correspond to the unique mean 
values of indefintely long records of each component. In the Ames data, averaging reduces the 
standard deviation of the random errors remaining in the measured values. 

The reference grid and a table indicating the number of measurements from reference 3 at 
each grid point is given in figure A2. At points with two measurements, averaging reduces the 
standard deviation of the independent errors, but retains the systematic errors present in both 
tables. At the remaining grid points in the common region one or both tables must be interpolated; 
this adds an approximation error to one or both values, and adds the average approximation error 
to the inaccuracy of the averaged measurements. This approximation error is small for sufficiently 
small interpolation intervals, or it may be small compared to the independent measurement errors. 
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For unfavorable conditions (large intervals and strong function nonlinearity) the averaging scheme 
may not provide the best accuracy at points with one measurement. 

The differences between the two sets of measurements are also given in figure A2 along 
with some statistics of the difference tables (means, standard deviations, and the table extremes) 
which describe the apparent random-measurement- error processes. All three force components 
and the yaw moment show small differences under 10Q lb, 30Q ft-lb, respectively. This excellent 
repeatability implies small random and interpolation errors. The roll- and pitching- moment 
differences are much larger; these exceed 25 Q ft-lb over a large part of the common domain and 
reach maximum values exceeding 100Q ft-lb. These two tables were also separated according 
to the number of measurements at each point (not shown) to verify that the largest differences 
occurred at grid points with two measurements and were not introduced by the interpolation. 
Further, these differences appear to be random with very small mean values; differences exceeding 
25 Q ft-lb are distributed throughout the grid, and no systematic behavior is discernible. 

The composite table is shown in figure A3. The reference grid for yaw has intervals of 10° 
to 15° for xp > 45°. This is sparse compared to the grid density below 45° and results in a coarse 
linearization of the aerodynamic functions with respect to xp in this region. Several of these 
(y, L, and Y M) are strongly nonlinear here, as can be seen in figure 2.3 of the text. This region 
is of particular interest in the case of single-point suspensions for which the MILVAN adopts a 
broadside attitude in equilibrium. The measurement grid density is greatest and the functions 
are best defined near xp = 0° which corresponds to the MILVAN’S equilibrium position with a 
two-point suspension. 


A2. Interference Tares 

Second, the interference tares must be determined and subtracted from the basic measure- 
ments. Reference 3 provides two tables of data obtained from yaw-sweep and pitch-sweep mea- 
suring sequences, respectively. These are quoted in figure A4 for all points in R 0 . 

The following factors affect the treatment. As seen from the grids in figure A4, the mea- 
surements do not cover all of R 0 ; the yaw-sweep data covers Rq for a < 20° and the pitch-sweep 
data covers R 0 for xp < 45°. In addition, the pitch-sweep data were taken without an image 
of the pitch-angle control arm so that these data measure the interference effects incompletely. 
Further, both tables were considered in reference 3 to be sufficiently erratic that their use was 
open to individual interpretation and was, therefore, left unsettled. The erratic measurements 
were attributed to random variations in the separation characteristics of the flow around the box. 
A review of the size of the yaw-sweep data indicates generally small values for drag and side-force 
interference, moderate values for lift and yaw moment, and large effects in roll and pitch moment. 
Thus, the interference tares are not small enough to discard as negligible effects. 

The first approach was the straightforward one of using the best available measurement or 
estimate; that is, the interference function was taken to be the yaw-sweep data for a < 20°, 
the pitch-sweep data for (20° < a < 45°, xp < 45°), and zero in the remainder of R 0 where 
no measurements were made. However, this approach 1) amplified the randomness of the basic 
measurements by adding the random measurement errors of the tare data, and 2) disturbed the 
systematic trends in the basic data by adding a function with distinctly different trends in the 
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three regions of Z 0 where it is differently defined. If the tare data are accurate measurements of an 
interference function which varies rapidly with (ip, a), then this treatment would tend to smooth 
the basic data. However, the results appeared to be a sum of independent random processes and 
are consistent with the hypothesis that the rapid variations in both sets of measurements are due 
to random measurement errors. 

In the approach adopted, the tare data were viewed as the superposition of a function which 
varies smoothly in (ip, a) and a random, independent measurement error; 

Af(ip,a) = Af(ip,oc) + rj (A3) 

where / refers to any aerodynamic component, and A/, A/ and r\ are the data, the actual inter- 
ference plus systematic measurement errors similar to those in the basic data, and the random- 
measurement errors with statistics similar to those of the random errors in the basic data, respec- 
tively. The data were therefore plotted and reviewed for systematic trends; that is, for consistent 
patterns versus ip, a with only a few extremes. 

The yaw-sweep data are plotted in figure A5. For the force components, the data plotted 
versus ip (fig. A5a) show consistent trends which can be estimated by averaging over a, but the 
data plotted versus a (fig. A5b) show no consistent trend nor any progression of curve shapes for 
different values of ip. For most values of ip, these curves can be interpreted as random measurement 
errors about a constant value independent of a. Straight-line least squares fits to the data showed 
nearly zero slopes at most values of ip. Further, the pitch-sweep data over the larger range, 
0° < a < 45°, did not indicate the presence of any significant a-dependent interference effect 
which might have been missed in reviewing the yaw-sweep data over 0° < a < 20°. Thus, the 
systematic variation of these tares can be estimated as independent of a and given by averaging 
the yaw-sweep data over a. The results, after minor smoothing of local curve reversals, are given 
in figure A6. The drag and side-force tares are seen to be small, below 10Q lb everywhere, while 
the lift interference reaches somewhat larger values at large yaw angles. 

Figure A6 includes some statistical properties of the difference from the original tare mea- 
surements; these differences are viewed here as the apparent random tare measurement errors, 
and also indicate the degree to which the tare data have been modified by our interpretation. The 
mean errors are zero by construction, and the standard deviations and extreme errors indicate 
measurement errors for the force components which are a moderate fraction of the tare data and 
small in absolute magnitude. Assuming the same random measurement error processes effect all 
sets of measurements then the statistics given in figure A6 should be close to those previously 
calculated in figure A2 from the differences in repeated measurements of the aerodynamics. A 
comparison of the two sets of statistics shows good agreement for the force components. 

The roll moment interference tare appears to be predominantly samples of a random vari- 
able. There is no discernible trend with ip (fig. A5a). Individual curves show reversals to local 
extremes and sign reversals at most grid points, and the data, when averaged over a, show the 
same properties again. Further, the data plotted versus a (fig. A5b) also show no discernible 
trend. Consequently, the roll-moment interference was estimated as the mean of the yaw-sweep 
data, —6 Q ft-lb, independent of either ip or a. This is a nearly negligible value compared to the 
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table extremes (— 72.4,63.6Q ft-lb). Thus, the roll moment interference tare table is interpreted 
to be almost entirely samples of a random process or random measurement errors, which is statis- 
tically consistent with the differences in repeated roll moment measurements given in figure A2. 
Evidently, the systematic roll moment interference is much smaller than the random measurement 
errors. 

The pitching-moment interference data shows a consistent trend vs xp, and a consistent, 
increasing trend with a at xp = 0° and low yaw angles. This trend declines to a nearly fixed 
value at xp = 45° and the data become erratic above 45°. Thus, systematic variations in both 
xp, a are discerned. For a < 20° the yaw-sweep data are used after smoothing of some local curve 
reversals. For the region a > 20°, the pitch-sweep data were reviewed; when plotted versus a 
these showed variations which 1) paralleled the yaw sweep data but were offset by a fixed bias, 
and 2) were constant for a > 20°. Consequently, the interference tare was estimated as fixed for 
a > 20°. The results are given in figure A6, and show large tares through most of k 0 . The error 
statistics indicate a smaller, moderate-sized error process (the large extreme error, 116Q ft-lb, 
seen in figure A6 is an isolated value in the error table) . 

Last, the yaw moment tares versus xp (fig. A5a) are nearly independent of a for xp near 0° 
and above 60°. Elsewhere, there is some consistent variation with a (fig. A5b) but the spread is 
small enough in magnitude to treat by averaging over a. The results are given in figure A6, and 
show modest tare magnitudes, and a random-measurement error process of about the same size. 

A3. Reference Table 

The basic data with the estimated interference tares subtracted are given in the text (fig. 2.2) 
along with plots versus xp and contour plots (figs. 2.3, 2.4). 
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( a ) Pitch-sweep data ; force components . 

Figure Al .- Wind tunnel data ; force and moment measurements . 
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(b) Pitch-sweep data; Moment components. 


Figure Al.- Continued. 
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( c ) Yaw-sweep data . 
Figure Al .- Concluded . 
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(a) Number of measurements each grid point. 

Figure A2.- Differences between yaw- and pitch-sweep data. 
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(b) Force differences. 
Figure A2 .- Continued. 
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- 31.5 

21.2 

- 9.2 

- 40.8 

- 70.3 

87.7 

27.9 

21.5 

3.8 

12.0 

- 18.0 - 25.3 

- 40.7 

- 53.6 

- 61.1 - 55.8 

- 14.0 

8.1 

6.3 

- 23.2 

- 62.3 

58.4 

22.9 

17.2 

32.4 

16.0 

- 14.4 - 29.7 

- 26.0 

- 9.9 

- 34.5 - 35.2 

- 51.4 

- 8.5 

- 8.1 

21.9 

18.8 

46.7 

23.7 

33.4 

42.4 

20.0 

- 9.8 - 23.6 

- 21.4 

- 8.7 

- 38.8 - 23.7 

- 55.9 

- 21.5 

- 25.3 

12.5 

32.2 

42.3 

3.3 

49.7 

28.7 

25.0 

11.8 - 21.8 

- 67.8 

• 

00 

3.9 - 9.5 

- 37.7 

- 22.3 

- 8.5 

6.1 

21.4 

25.1 

- 2.3 

19.8 

19.2 

30.0 

- 16.4 - 22.1 

- 46.3 

• 34.4 

- 18.7 27.6 

7.3 

- 12.9 

- 34.7 

- 39.5 

- 25.0 

- 16.1 

- 8.6 

3.5 

9.1 


a 

0.0 

- 1.2 

- 0.5 

0.9 

- 1.3 

- 5.1 

- 5.0 

YM/Q 

- 0.6 0.5 

1.3 

- 1.9 

0.7 - 16.6 

- 1.9 

4.9 

- 1.7 

2.0 

0.1 

- 2.9 

- 8.2 

- 7.6 

- 6.1 

- 4.8 

- 2.2 

- 0.6 

2.2 

- 0.5 

2.5 - 15.4 

- 4.8 

0.8 

- 7.7 

4.0 

- 1.0 

8.3 

0.3 

- 5.4 

- 5.1 

- 3.8 

- 2.1 

1.2 

2.2 

0.8 

4.7 - 13.0 

- 5.0 

- 0.9 

- 7.7 

6.0 

- 0.8 

5.1 

- 13.8 

- 26.5 

- 13.9 

- 7.1 

5.8 

5.4 

- 0.7 

- 0.5 

- 2.6 - 23.6 

- 2.5 

- 3.2 

- 7.7 

8.0 

1.4 

- 1.0 

- 17.8 

- 18.2 

- 12.5 

- 0.4 

3.2 

3.5 

0.1 

5.8 

1.9 - 26.9 

- 1.5 

- 3.9 

- 4.8 

10.0 

1.9 

1.0 

0.9 

5.4 

5.9 

10.4 

5.6 

- 8.7 

- 6.1 

1.2 

3.8 - 21.1 

2.3 

- 0.9 

- 4.4 

12.0 

2.8 

2.1 

6.8 

13.9 

22.6 

4.9 

7.1 

- 10.9 

- 3.6 

2.5 

8.8 - 14.5 

2.5 

- 0.7 

- 3.8 

16.0 

3.2 

2.1 

5.5 

3.4 

1.5 

8.5 

2.6 

0.9 

5.3 

0.0 

4.6 - 11.5 

0.1 

3.3 

3.0 

20.0 

- 1.9 

12.4 

12.4 

12.3 

28.1 

15.1 

17.0 

2.0 

7.1 

2.2 

3.9 - 13.1 

- 4.2 

0.8 

1.0 

25.0 

- 3.7 

5.9 

16.9 

14.8 

24.1 

25.2 

22.7 

- 1.2 

0.4 

2.2 

6.2 - 11.7 

- 5.5 

- 3.1 

- 0.9 

30.0 

- 3.2 

- 2.2 

0.1 

2.9 

0.9 

- 1.0 

1.0 

1.9 

- 0.6 

1.6 

1.8 - 13.2 

- 2.3 

- 5.0 

- 3.1 
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16.0 
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45.0 50.0 

60.0 

75.0 
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TABLE STATISTICS 



MIN 

MAX 

MEAN 

STD DEV 

RM 

- 109.7 

101.1 

5.5 

39.1 

PM 

- 70.3 

136.4 

- 8.8 

29.5 

YM 

- 26.9 

28.1 

-.2 

8.5 


( c ) Moment differences . 
Figure A2 .- Concluded . 
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a 

0.0 

60.3 

62.3 

70.0 

80.9 

91.9 

100.6 

D/Q 

113.9 128.0 

139.2 

149.8 

160.3 

172.2 

198.9 

213.6 

225.2 

2.0 

63.3 

64.6 

72.1 

82.1 

93.0 

102.1 

115.7 

129.8 

141.6 

151.9 

161.6 

174.4 

199.7 

212.9 

225.2 

4.0 

66.3 

68.1 

75.3 

83.4 

94.6 

104.9 

117.2 

131.7 

143.4 

153.7 

163.1 

177.0 

201.3 

213.2 

226.7 

6.0 

69.2 

71.8 

79.0 

86.7 

96.8 

107.1 

120.7 

133.9 

145.8 

155.7 

165.8 

180.1 

201.8 

212.4 

226.4 

8.0 

72.9 

75.1 

82.2 

90.5 

99.7 

109.8 

123.6 

136.6 

148.6 

158.8 

170.5 

184.3 

202.2 

212.2 

225.4 

10.0 

76.6 

78.4 

85.6 

94.1 

104.0 

114.4 

128.5 

140.2 

151.9 

162.6 

174.5 

186.9 

202.9 

212.6 

225.7 

12.0 

80.5 

82.0 

88.9 

97.8 

108.3 

119.7 

133.4 

144.4 

156.2 

166.5 

177.6 

188.8 

203.6 

211.6 

225.6 

16.0 

88.1 

89.0 

96.0 

105.6 

115.8 

128.6 

143.4 

153.7 

162.9 

172.7 

182.4 

191.5 

204.6 

210.1 

225.3 

20.0 

97.0 

99.1 

105.6 

114.7 

125.9 

137.4 

150.6 

158.8 

169.2 

178.0 

186.3 

194.0 

204.8 

208.6 

225.7 

25.0 

109.4 

111.1 

118.2 

127.3 

136.9 

148.2 

159.6 

167.9 

177.8 

186.0 

193.7 

199.1 

205.9 

206.6 

225.0 

30.0 

121.3 

122.4 

128.9 

139.4 

145.7 

156.4 

166.8 

176.9 

186.4 

194.1 

200.8 

203.9 

206.7 

206.1 

225.8 

35.0 

129.0 

133.4 

140.9 

148.2 

153.3 

164.7 

174.8 

188.8 

196.4 

202.3 

208.2 

208.3 

208.4 

204.2 

224.9 

40.0 

144.2 

145.3 

151.6 

158.0 

163.5 

170.5 

183.7 

196.6 

204.6 

211.3 

218.0 

215.3 

209.9 

204.3 

223.2 

45.0 

159.6 

161.8 

165.6 

167.4 

170.1 

180.4 

193.2 

204.0 

212.7 

219.1 

225.4 

221.4 

213.4 

201.2 

224.0 

a 

0.0 

- 4.0 

11.2 

20.7 

27.7 

32.2 

41.9 

Y/Q 

57.5 68.4 

73.6 

79.8 

82.2 

84.3 

89.4 

78.2 

3.1 

2.0 

- 3.0 

12.3 

21.5 

27.3 

31.2 

41.0 

56.5 

66.0 

71.6 

76.7 

78.9 

82.3 

90.3 

78.3 

3.5 

4.0 

- 2.5 

12.7 

22.1 

31.2 

31.2 

41.4 

55.7 

65.5 

70.5 

74.3 

76.1 

80.5 

90.8 

79.2 

5.1 

6.0 

- 3.0 

13.5 

24.7 

31.4 

33.2 

43.5 

54.7 

63.9 

68.5 

73.0 

76.4 

82.1 

91.9 

80.0 

4.4 

8.0 

- 3.8 

13.9 

26.6 

32.2 

34.6 

45.0 

54.8 

62.5 

68.3 

72.5 

76.2 

83.0 

91.9 

81.0 

4.6 

10.0 

- 4.4 

13.9 

29.4 

34.7 

38.9 

48.0 

55.6 

62.5 

66.8 

72.0 

76.5 

82.9 

92.1 

81.8 

4.9 

12.0 

- 5.8 

14.3 

31.0 

36.6 

42.0 

49.7 

55.9 

61.5 

67.0 

72.2 

76.4 

82.9 

93.2 

82.7 

4.6 

16.0 

- 7.3 

14.1 

34.5 

45.5 

52.8 

54.9 

53.2 

58.7 

63.5 

70.4 

75.6 

81.4 

92.2 

84.5 

5.4 

20.0 

- 11.4 

15.5 

34.2 

45.8 

55.3 

55.7 

51.9 

54.8 

60.2 

66.8 

72.5 

78.0 

88.8 

85.2 

4.5 

25.0 

- 12.4 

14.0 

31.0 

44.6 

54.2 

58.8 

56.8 

53.1 

59.3 

63.5 

67.9 

73.6 

85.3 

83.4 

5.1 

30.0 

- 12.5 

12.1 

25.7 

38.4 

47.9 

54.4 

56.9 

52.6 

57.4 

59.5 

63.0 

68.5 

80.5 

83.4 

3.5 

35.0 

- 5.4 

10.1 

17.4 

26.3 

34.8 

46.0 

48.8 

50.9 

51.9 

55.1 

58.3 

63.3 

73.4 

83.0 

- 0.6 

40.0 

- 3.2 

4.7 

15.2 

18.3 

25.3 

34.7 

41.1 

45.2 

47.3 

49.2 

51.0 

56.6 

67.7 

81.2 

- 1.2 

45.0 

- 0.5 

- 0.2 

6.5 

6.9 

13.8 

22.1 

32.7 

37.9 

38.9 

40.7 

42.5 

48.3 

59.8 

78.0 

- 6.6 

a 

0.0 

4.1 

3.6 

6.0 

8.4 

10.3 

16.3 

L/Q 

24.1 32.7 

39.2 

42.0 

39.3 

27.9 

5.3 

- 16.6 

- 5.8 

2.0 

8.7 

8.5 

10.8 

15.2 

18.6 

24.6 

33.5 

42.5 

49.3 

50.3 

45.1 

29.9 

2.1 

- 16.5 

- 5.2 

4.0 

12.3 

12.2 

17.4 

21.8 

27.3 

33.3 

42.4 

51.9 

59.6 

57.9 

49.6 

31.4 

0.0 

- 17.1 

- 6.2 

6.0 

15.5 

15.1 

20.1 

26.2 

33.1 

39.8 

51.6 

60.3 

66.2 

61.8 

50.3 

32.9 

0.6 

- 17.3 

• 4.8 

8.0 

18.4 

17.8 

21.1 

29.3 

37.1 

44.5 

57.9 

68.1 

72.2 

64.6 

49.5 

33.4 

1.3 

- 17.5 

- 4.5 

10.0 

19.6 

20.0 

24.2 

31.4 

39.6 

49.2 

61.6 

70.1 

73.5 

65.0 

50.5 

35.3 

3.0 

- 18.1 

- 4.2 

12.0 

22.9 

22.3 

26.7 

33.5 

42.8 

54.2 

65.1 

71.9 

69.7 

63.3 

52.5 

37.8 

4.8 

- 18.9 

- 3.9 

16.0 

26.6 

27.9 

33.4 

41.3 

49.5 

62.5 

76.8 

73.1 

68.2 

61.3 

53.4 

40.8 

8.4 

- 19.0 

- 3.5 

20.0 

34.3 

34.3 

39.3 

47.3 

55.0 

66.9 

77.1 

74.6 

69.9 

63.6 

56.0 

43.4 

12.0 

- 16.9 

- 3.0 

25.0 

43.3 

42.1 

45.9 

53.4 

58.3 

64.7 

72.7 

76.4 

74.3 

68.7 

61.5 

48.5 

17.4 

- 12.5 

- 2.5 

30.0 

52.1 

49.3 

51.3 

59.2 

62.5 

66.1 

68.9 

76.8 

76.6 

72.2 

66.0 

52.7 

22.5 

- 9.2 

- 1.9 

35.0 

53.9 

55.8 

57.7 

60.1 

67.1 

69.1 

72.9 

76.2 

77.6 

73.3 

68.9 

54.7 

26.4 

- 5.7 

- 3.3 

40.0 

60.1 

58.8 

62.3 

66.8 

68.5 

69.3 

73.7 

76.0 

75.1 

72.4 

69.7 

56.8 

31.0 

- 0.6 

- 1.1 

45.0 

61.9 

62.9 

65.2 

67.3 

66.9 

69.1 

72.6 

73.2 

73.0 

70.5 

68.0 

56.7 

34.2 

3.0 

- 1.9 

V- 

0.0 

4.0 

8.0 

12.0 

16.0 

20.0 

25.0 

30.0 

35.0 

40.0 

45.0 

50.0 

60.0 

75.0 

90.0 


( a ) Force components . 


Figure A3 .- Combined yaw - and pitch-sweep data . 
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RM/Q 


0.0 

7.3 

- 9.0 

7.3 

- 11.8 

20.2 

00 

1 

- 22.7 - 17.5 - 32.8 - 88.6 - 128.4 - 161.6 - 189.2 - 113.3 - 85.6 

2.0 

- 17.9 

- 12.0 

32.8 

- 2.8 

15.9 

13.7 

- 22.5 17.0 - 28.4 - 90.5 - 147.7 - 158.9 - 203.5 - 115.0 - 96.4 

4.0 

- 10.2 

5.9 

- 14.9 

- 46.3 

2.5 

- 7.3 

- 6.3 - 14.6 - 55.5 - 109.2 - 168.9 - 156.8 - 215.6 - 111.0 - 110.3 

6.0 

4.4 

- 2.7 

- 49.0 

- 32.9 

7.8 

- 0.8 

- 11.0 - 12.1 - 41.2 - 141.7 - 218.5 - 206.8 - 244.7 - 124.8 - 120.3 

8.0 

16.6 

- 10.9 

- 48.0 

- 29.2 

20.7 

17.4 

28.2 - 26.8 - 85.1 - 166.9 - 195.0 - 199.8 - 248.8 - 144.9 - 112.1 

10.0 

25.0 

- 3.8 

- 61.9 

- 19.8 

- 8.9 

9.7 

20.1 - 32.5 - 72.1 - 151.8 - 186.4 - 194.1 - 244.1 - 139.6 - 123.1 

12.0 

34.3 

0.0 

- 48.1 

- 3.4 

16.3 

- 8.9 

6.0 - 25.2 - 93.8 - 151.1 - 172.0 - 198.7 - 281.5 - 141.2 - 107.6 

16.0 

42.2 

12.6 

- 37.8 

- 26.8 

- 29.5 

- 18.5 

- 12.9 - 11.8 - 40.6 - 102.6 - 145.3 - 182.7 - 276.3 - 137.1 - 118.7 

20.0 

73.2 

2.6 

- 29.7 

- 2.2 

- 37.1 

- 40.8 

- 23.7 - 22.0 - 20.2 - 74.3 - 128.1 - 162.8 - 252.7 - 146.0 - 103.0 

25.0 

61.7 

31.6 

4.3 

- 5.9 

- 31.7 

- 29.8 

- 23.7 - 17.0 - 24.5 - 74.7 - 148.1 - 177.8 - 259.9 - 116.4 - 144.1 

30.0 

71.0 

51.5 

8.4 

33.8 

- 56.7 

- 46.0 

- 42.5 - 42.8 - 53.6 - 80.4 - 153.9 - 174.6 - 240.7 - 110.8 - 109.5 

35.0 

75.2 

19.2 

48.7 

77.5 

9.4 

- 71.8 

- 10.8 - 58.4 - 64.5 - 126.9 - 189.4 - 192.6 - 199.1 - 73.2 - 66.8 

40.0 

74.6 

39.4 

40.3 

- 22.7 

- 12.8 

- 22.5 

- 68.1 - 118.6 - 78.2 - 138.3 - 198.3 - 206.8 - 223.8 - 67.9 - 64.6 

45.0 

78.3 

64.9 

57.0 

23.4 

- 15.4 

- 26.9 

- 35.0 - 26.7 - 70.7 - 114.5 - 158.4 - 175.0 - 208.2 - 72.0 - 35.4 


PM/Q 


0.0 

- 47.0 

- 26.9 

- 33.8 

- 24.4 

7.3 

0.8 

- 5.5 

- 13.0 

- 13.4 

- 22.3 - 30.6 - 81.3 - 242.8 - 253.3 - 221.8 

2.0 

- 87.6 

- 62.3 

- 44.2 

- 17.1 

17.6 

17.8 

16.7 

23.6 

10.2 

- 6.2 - 18.5 - 73.0 - 243.2 - 262.6 - 209.0 

4.0 

- 93.1 

- 84.6 

- 57.0 

- 21.6 

21.7 

41.9 

39.9 

62.0 

40.4 

18.6 4.7 - 57.7 - 244.8 - 241.1 - 218.9 

6.0 

- 97.8 

- 95.6 

- 60.6 

- 13.1 

9.7 

31.0 

44.3 

70.9 

61.9 

39.9 18.1 - 101.6 - 253.4 - 249.6 - 221.5 

8.0 

- 104.7 

- 104.8 

- 59.8 

- 25.8 

4.4 

9.1 

30.4 

54.6 

67.8 

45.9 17.0 - 155.4 - 262.8 - 234.9 - 208.6 

10.0 

- 88.2 

- 88.0 

- 51.2 

- 30.9 

3.4 

- 0.2 

25.6 

53.4 

49.2 

6.4 - 37.6 - 170.4 - 248.2 - 245.3 - 213.6 

12.0 

- 60.1 

- 55.5 

- 36.4 

- 32.1 

- 5.2 

6.9 

40.3 

28.2 

28.6 

- 23.0 - 70.0 - 175.7 - 248.9 - 241.2 - 208.9 

16.0 

- 12.3 

- 9.4 

10.1 

- 3.0 

24.0 

26.5 

33.4 

- 17.4 

- 35.0 

- 76.9 - 102.1 - 162.8 - 244.9 - 220.6 - 222.9 

20.0 

38.2 

32.4 

50.5 

28.0 

46.5 

41.6 

34.4 

- 9.2 

- 35.4 

- 56.7 - 69.0 - 134.2 - 234.9 - 205.6 - 230.2 

25.0 

109.2 

81.9 

75.0 

66.2 

98.6 

60.3 

47.6 

8.8 

- 16.7 

- 28.9 - 41.3 - 104.1 - 212.4 - 211.9 - 235.4 

30.0 

155.3 

130.3 

133.5 

91.3 

118.1 

90.5 

74.3 

31.9 

- 19.6 

- 20.7 - 31.5 - 86.2 - 190.3 - 211.4 - 240.9 

35.0 

192.1 

184.3 

179.3 

149.4 

163.2 

111.1 

119.0 

49.5 

- 2.7 

- 4.9 - 7.1 - 56.3 - 154.6 - 176.0 - 263.1 

40.0 

202.9 

220.5 

162.4 

166.9 

151.7 

146.9 

127.2 

51.2 

50.5 

16.6 - 17.3 - 45.5 - 101.9 - 230.0 - 250.5 

45.0 

206.5 

190.3 

201.1 

188.5 

179.2 

150.1 

98.0 

63.2 

43.8 

20.0 - 3.8 - 34.6 - 96.3 - 130.3 - 176.0 


a 

0.0 

30.1 

- 60.5 

- 97.2 

- 68.8 

- 15.9 

38.4 

YM/Q 

111.3 166.1 

206.7 

241.7 

273.9 

287.7 

290.8 

239.4 

37.7 

2.0 

27.7 

- 52.6 

- 96.9 

- 70.5 

- 19.9 

34.6 

103.2 

156.4 

195.8 

231.1 

263.6 

280.8 

292.0 

238.6 

37.8 

4.0 

24.0 

- 37.9 

- 87.8 

- 68.1 

- 22.9 

31.1 

95.9 

148.1 

184.5 

220.3 

253.5 

274.4 

294.6 

238.9 

40.9 

6.0 

19.3 

- 15.3 

- 48.1 

- 45.4 

- 18.0 

28.0 

93.0 

141.0 

179.0 

212.2 

246.6 

273.5 

295.8 

238.4 

39.7 

8.0 

15.5 

5.8 

- 3.2 

- 8.0 

- 8.1 

29.8 

85.0 

130.8 

175.3 

208.0 

245.6 

276.3 

296.1 

238.7 

39.9 

10.0 

13.4 

15.0 

19.2 

22.5 

25.2 

46.9 

95.6 

134.4 

179.8 

213.7 

249.9 

278.0 

297.7 

238.3 

39.1 

12.0 

11.4 

23.7 

35.2 

45.5 

57.7 

55.8 

105.8 

143.0 

188.7 

222.3 

255.7 

280.2 

297.4 

236.6 

38.4 

16.0 

6.8 

40.0 

60.7 

77.7 

95.4 

80.8 

122.4 

168.3 
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Figure A3 .- Concluded . 
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Figure A4.- Interference tares. 


55 



a 


AD/Q 





ARM/Q 


0.0 

3.2 

2.0 

2.7 

8.4 

- 34.3 

5.2 

- 44.8 

61.0 

4.0 

3.2 

1.7 

4.5 

6.5 

- 7.7 

31.1 

24.3 

13.5 

8.0 

4.4 

1.7 

3.7 

6.5 

16.0 

- 9.8 

- 24.6 

- 32.9 

12.0 

2.6 

2.3 

3.7 

7.0 

- 44.7 

- 2.5 

- 21.5 

- 32.6 

16.0 

0.6 

2.1 

5.2 

8.9 

- 38.7 

54.7 

11.5 

- 27.0 

20.0 

1.7 

1.1 

4.4 

9 . 

- 28.1 

21.8 

60.0 

19.3 

25.0 

0.7 

- 0.8 

5.5 

5.5 

19.7 

- 3.1 

- 29.0 

- 25.1 

30.0 

0.3 

- 3.8 

1.5 

1.5 

- 8.4 

- 16.0 

- 62.6 

16.6 

35.0 

- 1.1 

- 6.7 

- 0.4 

5.8 

- 4.2 

- 41.2 

- 65.7 

- 13.9 

40.0 

- 4.3 

- 8.5 

- 7.4 

0.0 

29.3 

- 56.3 

- 44.9 

- 38.8 

45.0 

- 2.2 

- 3.4 

- 6.1 

- 2.6 

10.9 

- 95.9 

- 4.8 

- 35.5 

a 


AY/Q 




APM/Q 


0.0 

- 2.1 

4.1 

8.4 

11.1 

CO 

to 

1 

35.2 

75.8 

43.6 

4.0 

- 0.8 

1.9 

6.8 

7.7 

22.4 

53.4 

93.3 

41.1 

8.0 

- 1.9 

0.7 

1.2 

2.7 

45.9 

89.0 

70.4 

29.8 

12.0 

- 2.6 

3.4 

2.3 

4.4 

61.9 

55.5 

67.8 

75.0 

16.0 

- 0.9 

4.9 

3.9 

3.8 

57.1 

67.5 

74.7 

60.1 

20.0 

- 2.1 

5.1 

3.6 

4.7 

70.6 

73.5 

59.7 

80.4 

25.0 

- 1.1 

1.7 

1.8 

3.1 

67.7 

78.0 

115.9 

40.7 

30.0 

- 1.7 

- 2.3 

0.0 

0.6 

82.1 

62.1 

118.2 

51.3 

35.0 

- 1.5 

- 14.0 

2.6 

2.5 

89.3 

89.6 

72.5 

11.9 

40.0 

1.7 

- 14.3 

- 1.7 

5.1 

43.0 

- 12.0 

39.3 

34.8 

45.0 

1.0 

- 9.2 

- 0.6 

3.7 

75.1 

80.8 

49.6 

90.5 

a 


AL/Q 




AYM/Q 


0.0 

0.6 

2.5 

8.8 

6.4 

4.1 

40.0 

40.3 

27.2 

4.0 

- 0.2 

2.9 

2.8 

6.8 

0.3 

19.0 

28.2 

28.4 

8.0 

- 1.1 

1.2 

5.4 

8.6 

1.6 

1.5 

3.3 

- 9.5 

12.0 

- 3.2 

- 2.7 

0.4 

5.2 

0.4 

27.2 

11.0 

1.4 

16.0 

- 1.4 

- 3.5 

- 1.6 

2.5 

- 2.6 

26.5 

11.9 

13.9 

20.0 

- 2.0 

CO 

1 

- 2.3 

2.3 

- 2.5 

28.7 

11.3 

10.6 

25.0 

- 4.9 

- 8.0 

- 1.9 

2.4 

- 2.1 

- 4.0 

7.1 

2.8 

30.0 

- 4.4 

- 6.1 

- 6.0 

1.3 

0.2 

- 14.8 

1.6 

- 2.5 

35.0 

- 5.1 

- 1.0 

- 3.4 

0.5 

- 1.6 

- 29.1 

10.0 

- 2.2 

40.0 

- 5.5 

- 4.5 

<q 

CO 

i 

2.2 

0.7 

- 20.0 

1.1 

- 6.0 

45.0 

- 1.9 

- 0.1 

- 4.5 

- 3.2 

- 1.7 

- 5.1 

1.9 

- 8.1 

0 

0.0 

16.0 

30.0 

45.0 

0.0 

16.0 

30.0 

45.0 


( b ) Pitch-sweep data . 
Figure A4 .- Concluded . 
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APPENDIX B. -UNIVERSITY OF MARYLAND WIND TUNNEL DATA 

The measurements made at the University of Maryland (UOM) 8- by 11-ft wind tunnel are 
given in reference 6 for a smooth-surfaced 8- by 8- by 20-ft container, along with additional 
measurements over a smaller grid to test the effects of variations in box details (rounded ends and 
corrugated sides) and box size. 

Data from reference 6 for the smooth box are quoted in figure Bl. These include corrections 
for interference tares; the tares were measured only sparsely and at a < 25°, and, as in the ARC 
test, were reported to be erratic. The measurement grid in figure Bl is dense at low sideslip 
angle but empty in the interval 50° < /? < 90°. This provides good function definition at low 
headings in the vicinity of the trim attitude for the container carried with a two point suspension, 
but poor definition in the region where the side force and yaw moment have extreme magnitudes 
and strong nonlinearities, and in the vicinity of its broadside trim attitude when carried with a 
single-point suspension. In general, grid point density is a factor in the accuracy of estimating a 
nonlinear function from measurements. 

The following discussion considers; 1) the agreement of these data with the zero lines and 
mutual relationships given in section 3 from the symmetry of the box, and 2) the agreement of 
these data with the estimated static aerodynamics derived in section 4 from the Ames data. 

Note that the UOM data is taken at positive sideslip (negative ip). Therefore, equations 3.10 
and 3.12 of the text are combined to give a relation analogous to equation 3.12 for positive sideslip: 


where 


D \ 


D \ 

Y 


— cos <p w Y — sin <p w L 

L 


— sin <p w Y + cos <p w L 

RM 


-RM 

PM 


cos <f> w PM — sin <p w Y M 

\YMJ 

(/?,«) 

— sin <p w PM — cos <p w Y M J 


tan ai = tan /?/ cos a 
sin/?i = sin a cos/? 


(Bl) 


cot <f> w = tan a sin /I = tan ai sin/?i 

Further, throughout this appendix, Si,S2,— >$6> 8X6 understood to denote regions of Pi which are 
defined by changing ip to (3 in figure 3.3, and the Ames data at (/?, a) are understood to be obtained 
from the data given at (— /?, a) in figure 4.5 using the symmetry properties in equation 3.10. 


Bl. Drag 

First, the drag data are tested for the agreement of values at equivalent points in the regions 
Si and S 2 in accordance with equation (Bl) by calculating the discrepancy: 

6D{/3, a) = DU(p, a) - DU{(3 u <*i) (/ 3 , a) € Si or S 2 [B 2) 
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The letter “U” is appended in the nomenclature ( DU,LU,...YMU ) throughout this section to 
indicate the UOM data in figure Bl. The results, (fig. B3) show good agreement with discrepancies 
below 10Q lb throughout Si and S 2 . 

Second, these differences are eliminated by averaging the data at related points in Si and S 2 . 
The resulting revised drag table (omitted for brevity) can be compared with the Ames data by 
calculating the difference 

AD = D- DU (B3) 

where D,DU refer to the revised Ames data (fig. 4.5) and the revised UOM data, respectively. 
The superscript ( ) refers to revised data throughout this section. This comparison is included in 
figure B3 and shows small differences below 5% of the drag at all points in the UOM measurement 
grid. 


B2. Side Force and Lift 

The side-force data of figure Bl show good numerical agreement with its theoretical zero 
values at P = 0,90°; the discrepancies have the same sign at both extremes of P and this suggests 
a small tunnel bias in measuring side force. The lift data in figure Bl show good agreement with 
the expected zero value along the boundary a = 0 (unlike the Ames data) as well as at P = 90°. 
Only small revisions, using equations 4.9, 4.10 from the text, are needed to null these data tables 
along the expected zero-force lines. 

The revised data are omitted for brevity. Their agreement with equation Bl was tested by 
calculating the discrepancy from lift predicted from measurements at the equivalent point; 

6L(P,a) = LU((3, a) + sin <f> w YU{Pi t ai) - cos <f> w LU (Pi, oti) (P,a) e $1 or $ 2 (B4) 

The results are tabulated and plotted in figure B4; discrepancies are below 10Q lb over about half 
the region of comparison, but larger discrepancies, reaching 25Q lb, occur otherwise. The two 
estimates of lift in equation B4 are compared in the plot and this shows good to fair agreement in 
curve shapes. Further, 6L > 0 at all gridpoints; that is, the lift data is greater than lift predicted 
from measurements at the equivalent point. Lift can also be compared with its value predicted 
from the Y measurements alone; 

6L'(P,a) = LU(P,a ) + (YU(P, a) cot <j> w + YU(Pi,on)csc<f> w ) ( P,a ) G Si or S 2 (B 5) 

The results (not shown) were similar to those for 6L ; that is, SV exceeds 10Q lb over the same 
region as 6L above, and LU is greater than predicted from the side force data at all grid points 
in Si and 52- These differences are assumed to be due to errors in measuring lift and side force 
and these errors are at least as large as the discrepancies. 

A comparison with the estimate from the Ames data (fig. 4.5) is included in figure B4, where; 

) (Be) 

A L=L-LU J v ' 

These differences contain the estimation errors and the small effects of differences in wind tunnel 
model surfaces. The side-force differences exceed 10Q lbs at only a few points and are otherwise 
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small, and plots of Y and YU in figure B4 demonstrate this good agreement graphically. The lift 
results show significant differences in the range of 10 to 30Q lb over a large part of Ro- This result 
for A L is close to the results for 6L and 6L' just given above, both in magnitude of the difference 
and the region of significant differences. Further, A L is negative at all grid points so that LU 
exceeds the estimate from the Ames data throughout Ro. Collectively, these comparisons suggest 
that the UOM lift data is less accurate than the side-force measurements in the region of largest 
differences, with these differences indicating the probable size of the lift measurement errors. 

B3. Roll Moment 

A comparison of figures B2 and 2.3 shows a considerable difference from the Ames roll moment 
measurements. RM was assumed an antisymmetric function in section 3 on the basis of the UOM 
data. An antisymmetric function will satisfy 

RM((3, a) = -RM{p x ,a x ) 

RM = 0° at P = 0°, or (3 = 90° , or a — 0° or tan /3 = sin a 

The data (fig. Bl) reach table extremes of —125 Q ft-lb, 77 Q ft-lb. Although their magnitudes 
differ significantly, these extremes have opposite signs and occur at nearly equivalent points, 
consistent with equation B7. A symmetric function would have the same sign at equivalent 
points. In addition, the UOM data show small to moderate values along the boundaries where 
zero is predicted for an antisymmetric function in equation B7. These values are within the range 
of the discrepancies from the expected zero lines seen in the PM and YM data discussed next, 
and are assumed here to be measurement errors. The evidence does not rigorously exclude a 
small symmetric aerodynamic effect superposed with a dominant antisymmetric effect, but does 
indicate that such an effect, if present, is smaller than the measurement errors and would be 
negligible in its effects on load dynamics. 

As was done in section 4 , the data are separated into symmetric and antisymmetric compo- 
nents over the regions Si and S 2 in accordance with equation Bl, using: 

RMUS{P,a) =0.5{RMU{P,a) + RMU{p u a x )) \ , , 

RMUA{p, a) = 0.5 {RMU{P, a) - RMU{P x ,a x )) J 1 ’ 

Assuming roll is antisymmetric, then RMUS is the average measurement error and RMU A is the 
roll moment plus the antisymmetric component of the measurement errors. In the results (fig. B5), 
RMU A shows good agreement with the theoretical null lines. The average error, RMUS , exceeds 
RMU A in magnitude over much of Si and S 2 ; this occurs in the region where both RMUS and 
RMU A have small magnitudes (below 25Q ft-lb) and for this region RMU A is a very uncertain 
estimate of a small roll moment. Elsewhere, RMU S reaches 36Q ft-lb, but this is only a moderate 
fraction of RMU A in this region (50 to 90Q ft-lb). The data can be revised to agree with the 
theoretical zero lines and equation Bl. This revision provides an estimate of the roll moment 
which can be used in place of the null value adopted in the text. A plot of this estimate is 
included in figure B5. However, the estimated roll moment is small, below 50Q ft-lb at nearly all 
grid points and would have little effect on the static or dynamic behavior of the load compared 
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to the restoring torques applied by a suspension. Consequently, the approximation, RM = 0, is 
expected to be adequate for a realistic simulation of the suspended load. 

B4. Yaw and Pitching Moment 

The yaw-moment data in figures Bl and B2 show small to moderate nonzero values along 
the expected zero-yaw moment lines at /? = 0,90°. A comparison of figure B2 and 2.3 shows good 
qualitative agreement of curve shapes with the Ames data (after accounting for the sign reversal for 
yaw moment). The pitching-moment data also shows small to moderate nonzero values along its 
expected zero-moment lines at a = 0,/? = 90°. A comparison with figure 2.3 shows a considerable 
difference from the Ames data both qualitatively for (3 above 30° and quantitatively everywhere. 

The data can be revised to null the tables at the theoretical zero-moment lines using equa- 
tions 4.9, 4.10 (and replacing YC and LC in those equations with YMU,PMU, respectively), 
and then the discrepancies from values computed from data at the equivalent point are 

6PM = PMU - PMU\ _ 

PMU X (/?, a) = cos <t> w PMU{^ u a x ) - sm^YMUtfucn) 

The results for 6PM (fig. B6) show moderate discrepancies (up to 50Q ft-lb) over most of Si 
and S 2 with a small region of larger discrepancies. The pitching moment functions compared in 
equation B9 are also plotted separately in figure B6 and this comparison shows good agreement 
of curve shapes. 

The revised UOM and Ames data (fig. 4.5) are compared in figure B7; 

A PM = PM - PMU \ 

AYM = YM - YMU J [ ’ 

The numerical results show the same level of agreement for both components, with differences 
below lOOQ ft-lb at nearly all grid points, and below 50Q ft-lb over most of Ro- Plots are included 
and show good agreement of curve shapes for the two sources of data. 

Conclusions 

The UOM data show good agreement with all zero-force and zero-moment lines predicted 
by the analysis of section 3, and are consistent with the estimate of the aerodynamics given in 
section 4, particularly with the rejection of the Ames RM and PM data as being dominated by 
errors, and with the estimate of PM from the Ames yaw-moment data. 

Agreement with equation Bl can be tested in the regions Si and S 2 . Agreement is good in 
the case of the drag data. The lift and side-force data, and the pitch- and yaw-moment data 
show good agreement of curve shapes, but moderate to significant numerical discrepancies from 
equation Bl. 

The roll-moment data are consistent with an antisymmetric function, but the data indicates 
that roll is a small moment which can be neglected and that an estimate of RM from the data 
would have uncertainties of the size of the roll moment itself over most of the region Ro. 

Comparisons in the region Ro with the estimates derived from the Ames data show small 
differences in drag and side force, while lift, pitch, and yaw moments show good agreement in 
curve shapes, but moderate to significant numerical differences in some regions. 
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The UOM data lacks measurements in the interval 50° < < 90° so that nonlinear aerody- 

namics are poorly defined in this region and an extrapolation of this data to the region 54 using 
equation Bl would be similarly affected. 
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Figure Bl .- Wind tunnel data for a smooth-surfaced 8 - by 8 - by 20-ft container ; 
University of Maryland ( ref . 6 ). 
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(a) 6D/Q : Difference between drag measurements at equivalent points. 
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(b) A D/Q: Difference from revised Ames data. 

Figure B3.- Drag comparisons. Notes: 8D{P i a) = DXJ(p^a) - DU[pi t oti)] )9i(/?,a), a\{p } a) from 
equation Bl; and A D = D — DU . 
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(a) 6L/Q : Lift difference from lift derived from measurements at equivalent points. 
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(b) A Y/Q: Difference from revised Ames data. 
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(c) A L/Q: difference from revised Ames data. 

Figure B4.-Lift and side force comparisons. Notes: 6L = JaJ — LU 1 A ; 

LI 7 i (/ J , a ) = — sin (ft i ati ) +cos^ ti ,L£7(^i,a 1 ); and a ), ai (/ 0 , a ) from equation Bl. 
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Figure B4.- Concluded. 
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(b) RMA/Q: Antisymmetric component of roll moment data. 



& 

(c) RMU : Revised antisymmetric component 

Figure B5.- Roll moment decomposition. Notes: RMS = 0.5(RMU + RMUi)\ 

RMA = 0.5(RMU — RMUi)\ RMUi(/3,a) = RMU(/3i t ax) Pi{/3,ot) 9 ai(j9,a) from equation Bl. 
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(a) 6PM/Q: Difference from pitching moment derived from measurements at equivalent points. 
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(b) Comparison of pitching moment from measurements at equivalent points. 

Figure B6.- Comparisonof pitch and yaw moment data at equivalent points. Notes: 6PM = PMV—PMU\\ 
P MU i{f), a) = coa (fr^P MU (fa, ai) - nn4> w YMtJ oli)\ and <j> w {P,a), Pi(P,a), Qi(/9,a) from 
equation Bl. 
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( a ) Differences from revised Ames data . 

Figure B7 .- Pitch and yaw moment comparisons with Ames data . Notes : A PM = PM — PMU; 
and A YM — YM - YMU. 
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(b) PITCHING MOMENT 


(c) YAW MOMENT 



Figure B7.- Concluded. 
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APPENDIX C.-NORTHROP CORP. WIND TUNNEL DATA 

The measurements made at the Northrop Corp. 8- by 10-ft wind tunnel are given graphically 
in reference 4 for a smooth-surfaced box. Reference 5 contains limited additional data plus 
extensive data on a truck and tracked vehicle. The original tabulated data are apparently lost. 
The accuracy of the tables presented here in reproducing those data is limited by; 1) the resolution 
accuracy of the graphs; and 2) apparent revisions of the measurements in the graphical results to 
impose null values at some of the expected zero lines. 

The available data are the “lateral-directional” components (side force, roll, and yaw mo- 
ments) in both wind and body axes. Values of drag and pitching moment can be derived from 
these using the transformation relations between body and wind axes components 


Tb,w = E2(<x)Es(— 0) 


from which we obtain: 

D = (Y cos/? - YB)/ sm0 (Cl) 

PM = (RM cos a cos/? — RMB — Y\M sin a) /cos a sin/? (C2) 

where YB and RMB refer to the body-axes side-force and roll-moment components. These 
equations are singular at 0 = 0° and ill-conditioned at low /?, up to 10°, because of plot resolution 
limits. However, some of the affected drag values can be supplied from drag data taken at 0 = 0° 
(ref. 5) , while the PM values at /? = 0° can be filled in from the Y M data at a = 0° in accordance 
with the equivalence of these measurements discussed in the text. Lift data is available only at 
0 = 0°, and is omitted. 

The results are given in figure Cl. The coarse a-grid permits only limited comparison with 
the theoretical symmetry properties derived in section 3. A cursory examination shows; 1) good 
agreement of drag values along the boundaries a = 0° and /? = 0°, 2) close agreement of side 
force with the expected zero lines at 0 = 0° and 90°, 3) close agreement of roll moment with the 
zero lines of an antisymmetric function at tan^ = sin a, a = 0°, 0 = 0°, 4) good agreement of 
pitching moment with its expected zero lines at a = 0° and 0 = 90°, and 5) null yaw moment at 
0 = 0°, but a noticeable offset from the expected zero crossing at 0 = 90°. 
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Figure Cl.- NORTHROP wind tunnel data (refs. 4 and 5): Smooth-surfaced 8- by 8- by 20-ft container 
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data. Nevertheless, an estimate can be defined for this region consistent with the measured or theoretical 
values along its boundaries and the theoretical equivalence of points within the region. Discrepancies 
between the NASA Ames Research Center’s data and theoretical symmetry properties range from small 
to gross, depending on the component, and are removed to apply the extrapolation relations. These dis- 
crepancies are assumed to be due to measurement errors, but their sources could not be identified. Data 
from independent wind tunnel studies are reviewed; these are less comprehensive than the NASA Ames 
Research Center’s data but show good to fair agreement with both the theory and the estimate given here. 
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